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ABSTRACT 

The  purpose  of  this  research  was  to  develop  a  new  code  to  analyze  the  data  from 
AFIT’s  positron  annihilation  of  radiation  spectrometer,  PsARS,  system.  The  PsARS  was 
used  to  analyze  positronium  formation  at  the  surface  of  a  capillary  tube  coated  with  gold. 
Using  energy  and  geometry  filters,  the  PsARS  code  determined  the  ratio  of  ortho- 
positronium  to  para-positronium.  An  electric  field  applied  across  the  capillary  tube  wall 
forced  positrons  to  move  towards  the  gold  surface  of  the  capillary  tube  and  prevent  them 
from  dropping  into  the  Ore  gap  before  they  reach  the  surface.  Inside  the  Ore  gap,  the 
positrons  are  more  likely  to  form  positronium.  With  the  addition  of  the  electric  field,  the 
positrons  gain  energy  to  overcome  the  loss  of  energy  from  scattering.  Once  at  the 
surface,  a  secondary  electric  field  applied  from  the  capillary  tube  surface  to  the  vacuum 
chamber  wall  polarized  the  positronium,  increasing  their  lifetimes.  The  capability  to 
force  the  positrons  to  gather  on  the  gold  surface  of  the  capillary  tube  was  demonstrated. 
In  addition  to  the  PsARS  system,  a  positron  annihilation  lifetime  spectroscopy,  PALS, 
system  was  used  to  analyze  the  lifetimes  of  the  positrons. 
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POSITRON  ANNIHILATION  RATIO  SPECTROSCOPY  STUDY  OF 


ELECTRIC  FIELDS  APPLIED  TO  POSITRONIUM  AT  MATERIAL 

INTERFACES 

1.  Introduction 


1.1  Motivation 

The  study  of  positrons  is  important  to  the  DOD  because  they  are  valued  for  their 
ability  to  characterize  material  defects  without  destroying  the  material.  The  lifetimes  of 
positrons  injected  into  a  material  are  proportional  to  the  size  of  the  defects  in  the  material. 
This  allows  a  characterization  of  the  material  without  destroying  it  because  the  positrons 
annihilate  with  the  electrons  in  the  material,  leaving  the  material’s  atoms  or  molecules 
intact.  This  will  be  discussed  in  further  detail  in  chapter  two. 

Furthermore,  the  positron’s  use  as  a  possible  energy  source  in  space  travel  has 
lead  to  the  development  of  positron  confinement.  Researching  technology  that  is  used  to 
confine  positrons  is  important  because  large  volumes  of  positrons  could  be  built  up  to  be 
used  as  fuel  for  space  travel  or  to  create  a  gamma-ray  source.  If  antimatter  could  be  used 
as  a  practical  energy  source,  it  would  revolutionize  space  travel.  Per  nucleon,  nuclear 
fusion  releases  more  energy  than  fission,  but  no  process  has  been  developed  that  doesn’t 
consume  more  energy  than  it  produces.  After  atoms  undergo  nuclear  fission  or  fusion, 
they  leave  behind  daughter  products,  decreasing  the  efficiency  of  the  process.  However, 
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with  positrons,  when  they  annihilate  with  electrons,  100%  of  their  mass  is  converted  into 
energy.  Using  positrons  for  space  travel  would  increase  the  efficiency  of  the  fuel. 

The  limiting  factor  in  space  travel  is  the  consumption  and  storage  of  fuel  and  the 
fuel’s  energy  density.  For  a  mission  to  occur  within  a  reasonable  amount  of  time,  the 
velocity  of  the  craft  needs  to  be  high.  The  best  fuels  for  achieving  a  high  velocity  must 
have  a  large  energy  density  and  specific  impulse.  Antimatter  can  produce  energy 
densities  that  are  three  orders  of  magnitude  larger  than  fusion  or  fission  and  ten  orders  of 
magnitude  larger  than  chemical  reactions  [4].  In  an  anti-hydrogen  engine,  the 
annihilation  of  protons  and  anti-protons  release  pions;  the  pions  then  heat  a  propellant, 
increasing  its  specific  impulse  [6].  This  will  produce  a  specific  impulse  of  one  million 
seconds,  while  the  specific  impulse  of  a  common  chemical  propellant,  hydrazine,  is  only 
220  seconds  [6].  Besides  an  increase  of  specific  impulse,  antimatter  fuels  have  larger 
energy  densities,  too.  Chemical  fuels  have  specific  energies  of  on  the  order  of  10  J/kg 
whereas  the  antimatter  fuels  are  nine  orders  of  magnitude  larger  [18].  Specifically,  the 
energy  density  of  antimatter  is  9x10^^  J/kg  [18]  while  a  conventional  propellant,  such  as 

•y 

kerosene,  has  an  energy  density  of  3.5x10  J/L  [7].  Current  fuels  are  inefficient  and  thus 
the  majority  of  spacecraft  payloads  are  used  to  store  fuel  to  compensate  for  the  fuel’s 
inefficiency.  If  anti-matter  was  used,  depending  on  its  density,  it  would  reduce  either  the 
weight  or  amount  of  fuel,  allowing  more  room  for  equipment  and  for  the  craft  to  achieve 
higher  velocities,  potentially  decreasing  mission  travel  times. 

Besides  its  use  as  a  possible  fuel  for  space  travel,  past  research  on  antimatter 
propulsion  has  lead  to  the  invention  of  the  Penning  trap  that  confines  positrons.  A 
Penning  trap  confines  charged  particles  using  a  solenoidal  magnetic  field  to  confine  the 
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particles  in  the  radial  direction  and  static  electric  fields  to  confine  them  in  the 
longitudinal  direction  [4].  Once  in  the  trap,  the  particles  will  cool  by  releasing 
synchrotron  radiation.  In  an  experiment  conducted  at  CERN,  they  were  able  trap  one 
million  anti-protons  [4].  Furthermore,  a  variation  of  the  Penning  trap,  the  nested  trap,  can 
be  used  to  store  two  oppositely  charged  particle  groups  that  overlap  in  space.  This 
eliminates  the  need  of  having  separate  containers  for  both  particles,  saving  both  space 
and  money.  A  nested  trap  can  be  produced  from  five  cylindrical  electrodes  that  are 
biased  in  such  a  way  as  to  produce  an  elongated  potential  well,  see  Figure  1  [18]. 


+  -  + 


Figure  1.  Diagram  of  the  electrode  biasing  needed  to  produce  the  nested  trap  that  confines  both 
positrons  and  antiprotons.  Each  sequential  horizontal  electrode  is  biased  opposite  of  the  previous 
electrode. 

Research  on  positron  confinement  is  important  because  once  large  volumes  of 
trapped  positrons  are  created,  the  positrons  can  be  used  to  create  mono-energetic  gamma- 
rays  used  for  gamma-ray  modulation.  The  gamma-ray  laser  would  be  created  from  a  Bose 
Einstein  condensates  of  positronium  [23].  Current  spectroscopic  methods  for  studying 
different  materials  only  gather  information  from  using  ultraviolet,  visible,  and  infrared 
radiation.  The  only  methods  using  gamma-rays  come  from  positron  angular  correlation 
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annihilation  radiation  (ACAR)  and  Doppler  broadening  annihilation  radiation  (DEAR) 
measurements.  ACAR  measures  the  deviation  from  co-linearity  of  the  two  511  keV 
annihilation  photons  in  the  laboratory  frame  of  reference.  DEAR  measures  the  Doppler 
shift  of  the  two  511  keV  photons  [24] . 

The  gamma-rays  are  formed  from  the  annihilation  of  electrons  with  positrons, 
commonly  produced  from  radioactive  beta  decay.  Positron  confinement  would  reduce 
the  need  to  rely  solely  on  radioactive  decay  to  garner  gamma  rays.  Resides  using  electric 
and  magnetic  fields  to  confine  the  positrons,  those  fields  can  also  be  used  to  polarize  a 
quasi-stable  state  the  positron  forms  with  an  electron,  positronium,  increasing  its  lifetime. 

Resides  its  possible  future  use  as  an  energy  source,  the  study  of  positrons  is  also 
valuable  because  they  can  be  used  to  probe  solids  (metals,  ionic  solids,  and  molecular 
solids),  liquids,  and  gases  without  destroying  them.  Scientists  at  Eglin  AFE  are 
developing  a  positron  annihilation  radiation  spectroscopy  (PAS)  capability  to 
characterize  pore  distributions  in  explosives  [11].  Furthermore,  PAS  has  previously  been 
used  to  characterize  defects  in  solid  materials  and  semiconductors.  This  is  possible 
because  positronium,  or  Ps,  most  probabilistically  only  form  in  defects  or  gaps  in  solids 
with  a  periodic  lattice  structure  and  at  the  surface  of  materials.  Eased  on  the  information 
gained  from  measuring  the  annihilation  photons,  the  types,  size,  and  locations  of  defects 
can  be  derived. 

1.2  Previous  Student  Research  and  Objective  of  Current  Research 

The  objective  of  this  research  is  to  understand  how  thermal  Ps  interacts  with  different 
materials  and  at  different  interfaces  under  the  influence  of  an  inhomogeneous  electric 
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field.  How  the  lifetime  of  the  Ps  and  its  decay  mechanisms  are  influenced  by  external 
electric  fields  will  also  be  studied.  Ps  likely  forms  at  material  surfaces,  and  by  forcing  the 
Ps  to  form  at  the  surface,  its  interaction  with  the  electric  fields  will  be  studied.  This 
project  is  the  continuation  of  a  project  started  by  Dr.  Larry  Burggraf  and  Lt  Robert 
Slaughter  from  the  previous  year.  Lt  Slaughter  developed  a  positron  annihilation  ratio 
spectroscopy  (PsARS)  code  that  determines  the  ratio  to  which  Ps  will  annihilate  via  two 
or  three  photons.  The  second  part  of  Lt  Slaughter’s  work  involved  radioactive  Na-22 
deposited  inside  a  capillary  tube.  The  capillary  tube  was  chosen  for  its  ability  to  contain 
the  Na-22  and  moderate  the  positrons  as  they  travel  to  the  surface  of  the  tube.  The  goal 
of  this  research  is  to  discover  how  to  maximize  the  amount  of  Ps  formed  at  the  surface  of 
the  capillary  tube.  Once  at  the  surface,  Ps  annihilation  can  be  altered  through  the  use  of 
applied  electric  fields.  The  positrons  will  be  brought  to  the  surface  of  the  capillary  tube 
by  applying  an  electric  field  across  the  inner  and  outer  wall  of  the  capillary  tube  to  keep 
the  positrons  from  losing  too  much  energy  as  they  traverse  through  the  capillary  tube. 
Through  the  addition  of  an  external  electric  field,  the  Ps  may  become  polarized,  altering 
the  lifetime  and  decay  mechanism  of  the  Ps. 

Lt  Slaughter  developed  two  methods  to  collect  and  analyze  the  3yl2y  data,  the 
coincidence  and  non-coincidence  methods.  The  non-coincidence  method  collects  a 
spectrum  from  all  segments  from  the  detector  and  compares  the  energies  of  the  photons 
to  determine  the  number  of  two  and  three  photon  events.  The  counts  in  the  340  keV  to 
40  keV  region  are  compared  to  the  counts  in  the  511  keV  peak.  The  counts  are 
normalized  to  a  calibrated  standard  spectrum,  and  count  and  energy  factors  are  applied  to 
calculate  the  3yl2y  ratio.  In  the  coincidence  method,  developed  at  AFIT,  the  energy  and 
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timing  information  from  the  different  detector  segments  are  used  to  employ  a  geometry 
filter  to  distinguish  three  photon  counts  from  two  photon  counts.  Four  of  the  detector 
segments  look  for  three  photon  events  and  the  other  two  segments  look  for  two  photon 
events. 

Each  method  has  their  advantages  and  disadvantages.  The  advantage  of  the  non¬ 
coincidence  method  is  that  because  the  three  photon  counts  lie  in  the  Compton  continuum 
for  the  511  keV  photons,  the  background  is  reduced  when  the  number  of  three  photon 
counts  is  determined.  However,  with  no  information  on  which  segments  of  the  detector 
recorded  the  three  photon  events,  their  numbers  could  be  artificially  increased  from 
Compton  scattering  between  adjacent  segments.  The  advantage  of  the  coincidence 
method  is  that  with  the  geometry  requirements,  the  chance  that  three  photon  events  could 
result  from  Compton  scatters  is  reduced,  but  with  only  four  of  the  six  detector  segments 
looking  for  three  photon  events,  the  method  is  highly  inefficient.  By  combining  the  two 
methods  together,  a  new  method  will  be  created  that  combines  the  advantages  of  both  the 
energy  and  geometry  method.  Using  both  the  coincidence  and  non-coincidence 
methods,  Lt  Slaughter  measured  the  following  3yl2y  ratio  [2: 114]: 

Table  1.  PsARS  measurement  of  gold  coated  capillary  tube  at  1  atm  and  a  vacuum  level  of  less  than 

10'^  torr. 


Pressure 

Coincidence  Sy/ly 
(1/372) 

Non-coincidence 
3y/2y  (1/372) 

Vacuum 

2.44 

3.78 

STP 

2.92 

3.85 

These  PsARS  results  showed  that  Ps  was  forming  on  the  surface  of  the  capillary  tube. 
The  large  ratio  difference  between  the  two  methods  occurred  because  the  coincidence 
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method  requires  an  additional  geometrical  correction  factor  that  was  not  applied.  The 
absent  correct  factor  would  account  for  the  capillary  tubing  not  being  centered  inside  the 
vacuum  chamber.  The  non-coincidence  method  does  not  require  the  capture  of  all  three 
gamma  emissions  thereby  estimating  the  Sy/ly  ratio  better  [2:114].  Requiring  the  capture 
of  all  three  gamma-rays  is  not  required  because  the  method  uses  a  sample  with  a  known 
3y/2y  ratio  to  calibrate  the  method.  The  measurements  at  vacuum  are  lower  than  the 
measurements  at  STP  because  Ps  interacts  with  the  atoms  and  molecules  in  the  air  at  the 
surface  of  the  capillary  tube.  In  addition  to  the  3y/2y  ratio,  lifetime  measurements  of  the 
tube  confirm  that  Ps  is  forming  on  the  gold  surface  on  the  outside  of  the  tube. 

To  measure  the  lifetimes  of  the  Ps,  Lt  Slaughter  used  the  positron  annihilation 
lifetime  spectroscopy  (PALS)  method.  This  method  used  two  BaF2  detectors  to  measure 
the  gamma  rays  released  from  the  capillary  tube.  The  detectors  were  placed  under  the 
capillary  tube,  facing  each  other.  The  gamma-ray  released  when  the  Na-22  decays 
indicates  the  birth  of  the  positron.  When  the  positron  annihilates  with  an  electron,  two 
511  keV  gamma-rays  are  produced.  These  gamma-rays  indicate  the  loss  of  the  positron. 
The  difference  between  the  two  times  equals  the  lifetime  of  the  positron.  Once  a  lifetime 
spectrum  is  collected,  the  data  is  analyzed  using  a  program  call  PALSfit.  The  PALSfit 
program  uses  a  least  squares  regression  to  find  the  lifetimes.  This  method  will  be 
discussed  in  further  detail  in  chapter  3.  Figure  2  shows  the  lifetime  spectrum  collected 
from  measuring  the  capillary  tube  for  one  week  [2: 117]. 
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Lifetime  (ns) 

Figure  2.  Lifetime  spectrum  collected  from  the  gold  capillary  tube  using  the  digital  PALS  system. 

The  measurement  lastly  approximately  one  week  [2:117]. 

This  figure  shows  that  several  long  lived  Ps  species  exist,  probably  at  the  surface  of  the 
capillary  tube.  However,  these  peaks  may  also  be  the  result  of  noise.  Due  to  the 
significantly  low  number  of  counts,  the  PALS  fit  software  was  not  able  to  find  a  best  fit  of 
the  lifetime  spectrum.  To  estimate  the  intensity  of  the  predominate  Ps  term  (47  nsec),  the 
number  of  counts  in  the  47  nsec  term  is  divided  by  the  total  number  of  counts  and 
multiplied  by  their  relative  lifetimes.  Lt  Slaughter  calculated  the  intensity  of  the  47  nsec 
lifetime  to  be  approximately  8.5%  [2:117]. 

Although  many  objectives  were  met  by  Lt  Slaughter  with  his  development  of 
AFIT’s  PsARS  code,  numerous  improvements  can  be  made.  When  Lt  Slaughter  heated 
the  capillary  tubes  to  deposit  the  NaCl  solution,  much  of  the  sodium  leaked  out.  To 
decrease  the  amount  that  escapes,  a  possible  improvement  to  the  method  would  be  to  heat 
both  ends  of  the  tube,  at  a  lower  heat,  and  allow  for  the  solution  to  evaporate  over  a 


longer  period  of  time.  Also,  in  the  coincidence  method,  Lt  Slaughter  only  used  four 


channels  to  record  the  three  gamma  spectrum.  If  all  six  were  used,  the  measurement  time 
would  have  been  greatly  reduced.  Part  of  this  research  will  involve  configuring  all  six 
channels  to  communicate  together.  Lt  Slaughter  was  forced  to  only  use  four  channels  to 
look  for  three  photon  events  because  the  electronics  were  malfunctioning  due  to  an 
unknown  reason  that  prevented  all  six  from  communicating  together.  Furthermore,  he 
only  achieved  a  vacuum  pressure  of  less  than  10’  torr,  with  the  roughing  pump.  To  use 
the  ion  pump,  a  starting  pressure  of  less  than  7.5x10’"^  torr  is  needed.  With  the  ion  pump, 
pressures  less  than  10’  torr  will  be  achieved.  At  these  lower  pressures  there  are  fewer 
atoms  and  molecules  for  the  positrons  to  interact  with  at  the  surface  of  the  capillary  tube. 

Positron  research  is  important  to  the  DOD  because  they  can  characterize  materials 
without  destroying  the  material.  Furthermore,  the  size  distribution  of  the  defects  can  be 
calculated  from  analyzing  the  positron  annihilation  lifetime  spectrum.  In  addition, 
research  on  positron  confinement  and  lifetime  extension  will  lead  to  the  creation  of 
gamma-ray  lasers  from  Bose  Einstein  condensates  of  positronium.  These,  and  further 
details  of  positron  and  positronium  theory  will  be  discussed  in  the  next  section. 
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2.  Theory 


The  positron  is  the  antiparticle  of  the  electron.  It  has  the  same  mass,  spin,  and  charge 
as  the  electron,  but  is  positively  charged.  In  a  vacuum,  the  particle  is  stable  with  a 
lifetime  greater  than  4x10  years.  However,  when  a  positron  travels  through  matter,  it 
will  annihilate  having  a  lifetime  that  is  inversely  proportional  to  the  local  electron  density 
[1:4].  Upon  annihilation,  the  masses  of  the  electron  and  positron  are  converted  to  energy 
which  is  released  as  gamma  rays.  The  two  most  common  numbers  of  gamma  rays 
released  are  two  and  three.  Higher  numbers  of  gamma  rays  can  be  released,  but  are  very 
improbable.  Furthermore,  a  radiationless  process  of  electron  emission  and  a  single 
gamma  ray  emission  can  occur,  but  because  they  can  only  occur  in  the  presence  of  a 
nucleus  to  conserve  energy,  these  events  are  very  improbable.  [1:4-6]. 

The  most  probable  annihilation  occurs  when  the  spin  of  the  two  particles  are  anti¬ 
parallel,  resulting  in  the  singlet  spin  state.  This  process  releases  two  gamma  rays.  The 
cross  section  of  this  reaction  is 


+  Ay +  \  — 7 

= 7  [  2  ,  +  Vr  - 1 )  ■ 

y+\  r  -1 


r-1 


where  ro  is  the  classical  radius  of  the  electron  and  y  is  the  Lorentz  factor,  which  is 
calculated 


(2) 


where  P  = 


V 

c 


(3) 
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with  V  equaling  the  speed  of  the  positron  and  c  equaling  the  speed  of  light  in  a  vacuum 
[1:5].  In  most  situations,  the  positron  is  traveling  at  non-relativistic  speeds,  and  the  cross 
section  equation  reduces  to 


(4) 


The  rate  at  which  the  electron-positron  pair  decays  into  photons  is  referred  to  as  the 
annihilation  rate.  It  is  calculated: 


B,p  =  I  ^SM  P))5\r^p){<ys,M  p) 


(5) 

(6) 


where  as,M  is  the  annihilation  cross  section,  \\f  is  the  wave  function  of  the  electron  and 
positron  pair,  rep  is  the  average  distance  between  the  positron  and  electron,  and  Ds  is  a 
constant  based  on  the  spin  [3:19].  In  equation  5,  the  sum  is  over  all  electron-positron 
pairs  in  the  system  and  in  equation  6  the  sum  is  where  the  spin,  S,  equals  0  andl.  Each  Ds 
is  calculated 


«  7.989^5  ' 

(7) 

Dj«7.040//5 

When  the  total  spin  equals  zero  (S=0),  the  positron  electron  system  is  in  the 
singlet  state  and  when  the  spin  equals  one  (S=l),  the  system  is  in  the  triplet  state  [3:19]. 
If  the  spins  are  anti-parallel,  and  the  annihilation  occurs  when  the  positron  has  low 
energies,  two  annihilation  photons  are  released,  each  with  approximately  511  keV’s  of 


energy.  In  the  laboratory  reference  frame,  the  photons  are  released  almost  180  degrees 
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apart.  In  the  center  of  mass  frame,  the  photons  are  released  exactly  180  degrees  apart.  If 
the  spins  of  the  electron  and  positron  are  parallel,  three  gamma  rays  are  usually  released 
[1:5,7].  When  three  photons  are  released,  they  can  have  a  spectrum  of  energies  from  zero 
toSllkeV  [25]. 

Whereas  conservation  of  momentum  and  energy  force  two  photon  annihilations  to 
have  the  same  energy,  released  collinear,  the  extra  degree  of  freedom  in  the  three  photon 
annihilation  allow  the  photons  to  have  a  spectrum  of  energies  and  angular  distributions 
[25].  In  a  triple  coincidence  experiment  run  by  DeBenedetti  and  Siegel  at  the  Carnegie 
Institute  of  Technology,  the  energy  of  the  three  Sy  photons  were  measured  at  varying 
angles.  They  used  three  Nal(Tl)  scintillation  detectors,  2.5  cm  thick  and  4  cm  in  diameter 
placed  12  cm  from  a  Na-22  source  [27].  Figure  3  shows  the  energies  of  the  three  photons 
at  different  angles  from  the  source  [reproduced  from  reference  27] .  Also  included  is  the 
case  of  two  511  keV  photons  as  an  energy  reference  for  each  detector. 
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[27]. 


With  each  Nal(Tl)  detector  located  120  degrees  from  the  source,  each  photon  is  going  to 
have  an  energy  of  340  keV.  As  the  angle  between  the  detectors  increase,  as  shown  in  the 
second  case  of  150  degrees,  the  photon  takes  away  less  energy,  270  keV,  and  the  smallest 
angle,  60  degrees,  took  away  the  remaining  475  keV.  However,  every  photon  that  is 
released  at  120  degrees  does  not  have  340  keV.  As  shown  in  the  last  case,  each  photon 
was  released  at  a  different  angle,  and  all  three  had  differing  energies.  The  energy  of  a 
photon  is  dependent  upon  the  angle  at  which  the  other  two  photons  are  released.  Because 
there  are  essentially  infinite  combinations  of  angles  with  which  the  photons  can  be 
released,  each  photon  can  have  a  spectrum  of  energies  from  zero  to  5 1 1  keV. 

The  most  common  source  of  positrons  is  when  Na-22  undergoes  beta  plus  decay. 
In  beta  plus  decay,  a  proton  is  converted  into  a  neutron  and  a  positron  and  neutrino  are 


13 


emitted.  When  Na-22  decays,  90.4%  of  the  time  it  decays  via  to  an  excited  state  of 
Ne-22.  3.7  psec  after  the  decay,  a  1.274  MeV  gamma  ray  is  released.  Figure  4  details 
the  decay  scheme  and  equation  8  shows  the  decay  reaction  [12:7]. 


22. ,  2.602  a 

Na  - 


(f  (0.1  %) 


22 


Ne 


(90.4  %) 
EC  (9.5  %) 

y  (1 .274  MeV) 


Figure  4.  Decay  scheme  of  Na-22  [12:7]. 


f^Na  llNe  +  e^  +y  +  v  (8) 

where  y  is  the  gamma  ray  and  v  is  the  neutrino.  Due  to  the  fact  that  the  neutrino  can  be 
released  with  a  spectrum  of  energies,  the  positron  is  also  released  with  a  spectrum  of 
energies,  see  Figure  5  [2:6]. 


(KeV*) 

Figure  5.  Positron  energy  spectrum  from  Na-22  beta  decay  [2:7]. 
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Before  annihilation  occurs,  the  positron  and  electron  can  form  a  quasi-stable 
bound  state,  called  positronium  (Ps).  Ps  has  a  reduced  mass  that  is  half  of  the  hydrogen 
atoms  reduced  mass.  With  a  reduced  mass  that  is  half  of  hydrogen’s,  Ps  has  a  ground 
state  binding  energy  (6.8  eV)  equal  to  half  of  hydrogen’s  ground  state  binding  energy 
(13.6  eV).  If  the  singlet  state  is  formed,  the  Ps  is  called  para-positronium  (p-Ps),  while 
the  triplet  state  is  called  ortho-positronium  (o-Ps).  The  number  of  gamma  rays  released 
for  para  and  ortho  is  given  by 

(-1)"'^  =(-1)^+^  (9) 

where  S  is  the  spin  of  the  Ps,  L  is  the  orbital  angular  momentum  of  the  Ps,  and  Uy  is  the 
number  of  gamma  rays  released.  In  the  ground  state,  p-Ps  can  only  emit  even  numbers  of 
gamma  rays  while  o-Ps  can  only  emit  odd  numbers  of  gamma  rays.  P-Ps  mostly 
commonly  emits  two  gamma  rays  while  o-Ps  mostly  commonly  emits  three.  Like  free 
positron  annihilation  with  an  electron,  higher  numbers  of  gamma  rays  can  be  released, 
but  these  situations  are  very  unlikely.  The  rate  of  ortho  to  para  formation  is  3:1. 
Therefore,  three  gamma  annihilation  is  more  abundant  for  Ps  than  for  free  positrons  [1:7- 
8].  In  a  vacuum,  Ps  can  only  annihilate  via  three  photons  to  conserve  spin  angular 
momentum  [12:81]. 

The  cross  section  ratio  for  three  gamma  to  two  gamma  annihilation  is  1/372  [14]. 
The  ratio  is  calculated 


_  ^^2 
O" ^ 


(10) 


where  03  is  the  cross  section  of  the  o-Ps,  02  is  the  cross  section  of  the  p-Ps,  13  is  the 
lifetime  of  the  o-Ps,  and  X2  is  the  lifetime  of  the  p-Ps.  The  factor  of  three  is  necessary 
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because  o-Ps  is  three  times  more  likely  than  p-Ps  to  form  [26].  Given  that  the  intrinsic 
lifetime  of  p-Ps  is  0.125  ns  and  142  ns  for  o-Ps,  the  cross  section  ratio  is  0.00268,  very 
close  to  the  known  1/372  ratio  [12:81]. 

In  a  gas,  or  at  surfaces,  Ps  is  formed  when  a  hyper-thermal  positron  collides  with 
an  atom  or  molecule  and  takes  away  one  of  its  electrons, 

+X^Ps+X^  (11) 

where  X  is  the  gas  atom  or  molecule.  For  Ps  to  form. 


2 

np. 


where  Eps  is  the  energy  of  the  positron,  Ei  is  the  ionization  energy  of  the  atom  or 
molecule,  and  nps  is  the  principal  quantum  number,  all  numbers  are  in  eV  [1:27,151].  In 
a  solid,  Ps  can  form  in  defects  in  the  material,  in  the  natural  voids  in  the  material’s 
matrix,  or  at  material  interfaces  [1:27].  At  material  interfaces  where  the  lattices  do  not 
match  up,  defects  occur  that  act  as  positron  traps  [12:78].  Furthermore,  the  arrangement 
of  the  atoms  in  the  material  can  create  potential  wells  where  the  Ps  can  become  trapped 
[1:27]. 

Inside  pure,  single  crystal  metals,  Ps  cannot  form  and  persist  due  to  the  high 
density  of  the  electrons,  but  can  form  at  the  surface  of  the  metal  [8].  As  the  positron 
approaches  the  surface,  it  can  reflect  back  due  to  the  surface  potential  barrier,  escape  as  a 
positron,  escape  as  Ps,  or  become  trapped  at  a  surface  state.  If  the  positron  work  function 
is  positive  and  the  Ps  work  function  is  negative,  the  positron  will  capture  an  electron  and 
escape  as  Ps  [8].  The  Ps  will  have  a  kinetic  energy  less  than  or  equal  to  the  negative  value 
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of  the  Ps  formation  potential,  eps.  The  formation  potential,  in  eV’s,  can  be  calculated 
from: 


6.8 

=9^ +  9-  — r 

np,. 


(13) 


where  (p+  and  (p.  are  the  positron  and  electron  work  functions,  respectively  [1:27].  For 
ground  state  Ps,  in  most  metals,  the  formation  potential  is  negative,  allowing  for  Ps  to 
leave  the  surface.  For  most  cases,  the  Ps  will  be  in  the  ground  state,  nps=l,  and  the 
binding  energy  of  the  Ps  will  be  6.8  eV  [1:27]. 

If  the  work  function  of  the  positron  is  negative,  the  positron  will  accelerate  toward 
the  barrier  and  if  it  has  enough  energy,  it  can  escape  from  the  material.  However,  positron 
emission  is  suppressed  as  the  temperature  decreases,  with  no  emission  as  the  temperature 
reaches  zero  Kelvin  [8].  A  negative  positron  work  function  allows  for  spontaneous 
positron  emission  from  the  surface  [8].  A  negative  work  function  is  defined  as  the  energy 
of  the  positron  which  has  been  thermalized  as  it  travels  through  the  material  [9].  This 
means  that  the  positron  will  leave  the  material  with  a  minimum  amount  of  kinetic  energy 
[28]. 

Gold  is  one  material  with  a  negative  positron  work  function  [9].  In  this 
experiment,  because  the  capillary  tube  is  coated  with  gold,  it  is  expected  most  of  the 
positrons  ejected  from  the  surface  will  possess  an  energy  equal  to  the  positron  work 
function  of  gold,  approximately  negative  1  eV  [9] .  However,  this  only  applies  if  the  gold 
surface  is  clean.  In  this  experiment,  the  gold  surface  is  coated  with  multiple  layers  of 
hydrocarbons  and  moisture  that  build  up  from  exposure  to  air.  These  hydrocarbon  layers 
will  have  an  effect  on  Ps  formation  and  annihilation,  altering  the  experimental  results. 
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The  study  of  negative  work  function  materials  is  important  because  mono-energetic 
positron  beams  can  be  created  from  them.  Further  studies  can  be  conducted  by  coating 
the  gold  with  another  material  and  the  positrons  emitted  at  the  gold  boundary  layer  can  be 
used  to  investigate  the  defect  characteristics  of  the  added  material.  In  addition  to 
studying  interactions  in  another  material,  a  secondary  applied  potential  between  the 
vacuum  tube  and  capillary  tube  intends  to  be  used  to  suppress  the  positron  emission  or 
polarize  the  Ps  on  the  gold  surface  to  influence  lifetimes. 

Positrons  are  commonly  formed  through  beta  plus  decay,  with  energies  of  several 
hundred  keV.  Once  the  positron  is  formed,  it  will  lose  energy  through  scattering. 
Although  rare,  a  positron  can  lose  half  of  its  energy  in  a  single  head-on  collision  with  a 
molecule  or  atom  in  the  material  the  positron  is  traveling  through  [3:117-118].  Most 
scatters,  though,  are  glancing  collisions,  resulting  in  an  average  energy  loss  of  30-50  eV 
per  collision,  with  a  maximum  of  100  eV  [3:117-118].  The  positron  will  most  likely  form 
Ps  if  its  energy,  W,  lies  within  the  Ore  gap.  The  Ore  gap  is  defined  as  the  energy 
interval 

(14) 

or 

lG-^<W<Ia  (15) 

where  1g  is  the  first  ionization  potential  of  the  molecule,  R/2  is  the  binding  energy  of  Ps 
in  a  vacuum,  and  Wex  is  the  electronic  excitation  threshold  of  the  molecule.  When  the 
energy  of  the  positron  falls  below  the  Ore  gap,  it  does  not  have  enough  energy  to  ionize 
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the  molecule  and  free  an  electron  to  form  Ps.  If  the  energy  of  the  positron  exceeds  the 
Ore  gap  level,  ionizations  and  electronic  excitations  will  dominate,  but  Ps  formation  is 
still  possible  [3:123]. 

The  rate  at  which  Ps  annihilates  is  called  the  annihilation  rate  and  is  influenced  by 
the  defects  and  pores  inside  the  material,  and  through  quenching.  Quenching  is  the 
process  through  which  o-Ps  is  transformed  into  p-Ps  [3:3].  With  quenching,  the 
annihilation  lifetimes  decrease  because  the  longer-lived  o-Ps  lifetimes  are  replaced  with 
the  shorter- lived  p-Ps  lifetimes.  The  annihilation  rate  of  the  Ps  is  dependent  on  the  size 
of  the  trap.  Assuming  the  trap  is  spherical,  with  a  radius  R,  the  annihilation  rate  for  o-Ps 
is 
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where 


/?o(«m)  =  R  +  .166  (17) 

for  ionic  and  molecular  solids  [3:90-91].  In  polymers,  the  p-Ps  will  decay  quickly  when 
the  positron  annihilates  with  its  paired  electron,  also  known  as  self  annihilation.  The 
longer-lived  o-Ps  can  decay  with  its  paired  electron  or  via  the  pick  off  method,  where  the 
positron  decays  with  an  electron  from  a  neighboring  atom  or  molecule  [3:3].  In  addition, 
a  second  process  called  exchange  quenching  can  occur.  This  is  where  the  o-Ps  exchanges 
its  paired  electron  with  an  electron  from  a  nearby  atom  or  molecule,  and  transforms  into 
p-Ps,  where  it  will  quickly  decay.  A  third  process,  called  chemical  quenching,  occurs 
when  the  o-Ps  undergoes  a  chemical  conversion.  Two  such  examples  are: 
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o-Ps  +  X  —>0  -  PsX  +  energy 
o-Ps  +  XY  ^o-PsX+Y 


(18) 


where  X  and  Y  are  atoms  or  molecules  [1:328-329].  All  of  these  processes  will  shorten 
the  o-Ps  lifetime  from  142  ns  to  a  few  ns  [3:3].  Through  the  application  of  positron 
annihilation  lifetime  spectroscopy  (PALS)  on  a  material,  the  o-Ps  lifetime  can  be 
determined  and  using  equation  16,  the  size  of  the  defects  in  the  material  can  be 
characterized. 

PALS  can  be  used  to  measure  both  the  size  and  concentration  of  defects  in  a 
material.  The  longer  lifetime  components  of  a  PALS  spectrum  are  due  to  a  decrease  in 
the  local  electron  density.  The  strength  of  the  longer  lived  component,  or  intensity,  is 
related  to  the  concentration  of  the  defects  in  the  material  [19].  Lifetime  measurements 
are  commonly  conducted  by  measuring  a  positron  source,  Na-22.  The  source  is  placed 
directly  between  the  two  detectors  used  to  conduct  the  measurements.  Detectors  with  fast 
timing  components  are  needed  to  be  able  to  resolve  small  time  resolutions.  Because  the 
1.274  Mev  gamma  ray  from  the  Na-22  is  released  almost  simultaneously  with  the 
positron,  it  indicates  the  birth  of  the  positron.  The  time  for  the  positron  to  thermalize,  a 
few  picoseconds,  is  so  small  compared  to  the  lifetime  of  the  positron  that  it  can  be 
neglected  [19].  Once  thermalized,  the  positron  diffuses  through  the  material  before  it  is 
trapped.  The  diffusion  length  determines  how  many  atoms  the  positron  can  probe  for 
defects  before  it  annihilates  [19].  The  diffusion  length  is  thus  an  indication  of  how  well 
the  positron  can  probe  the  material  for  defects.  When  the  positron  annihilates,  the 
annihilation  photon  indicates  the  positron’s  death.  The  difference  between  the  two  times 
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at  which  the  birth  gamma  and  annihilation  photon  are  detected  equal  the  lifetime  of  the 


positron  [19]. 

When  using  PALS  to  determine  the  lifetimes  of  the  positrons,  the  source  activity 
must  be  low.  If  the  activity  is  too  high,  more  than  one  positron  will  be  present  in  the 
material  at  the  same  time  and  the  start  and  stop  signals  could  be  mixed  up.  However,  if 
the  source  activity  is  too  low,  the  measurement  will  take  longer  to  collect  [19].  To 
produce  a  complete  spectrum,  one  million  counts  are  needed. 

Annihilation  rates  can  be  altered  through  the  use  of  external  electric  or  magnetic 
fields.  Applying  an  external  magnetic  field  to  the  Ps  will  cause  the  magnetic  subset  of 
the  ortho  and  para-positronium  to  mix,  increasing  the  lifetime  of  the  p-Ps  and  decreasing 
the  lifetime  of  o-Ps  [3:91].  An  experiment  conducted  by  Stevens  and  Lichtenberger  at 
the  University  of  Guelph,  Ontario,  Canada,  showed  that  the  lifetimes  of  positrons 
annihilating  in  polyethylene  were  extended  through  the  use  of  electric  fields  [16].  The 
experimental  setup  consisted  of  the  Na-22  sandwiched  between  two  nickel  foils.  Each 
layer  on  the  foil  was  followed  by  a  layer  of  polyethylene  and  then  a  layer  of  brass. 
Varying  high  voltages  were  applied  to  the  brass  layers  to  produce  static  electric  field 
across  the  polyethylene  [16].  Table  2  shows  the  results  [reproduced  from  reference  16]. 


Table  2.  Experimental  values  of  the  average  lifetimes  and  intensities  of  the  positron  annihilation  in 
_ _ polyethylene  under  various  static  electric  field  values  116]. 


E 

(kV/cm) 

B1 

H 

Hf 

h 

(nsec) 

T2 

(nsec) 

■^3 

6isec) 

T4 

(nsec) 

0 

23.5 

46.6 

10.6 

19.3 

0.145 

0.352 

0.663 

2.39 

30 

19.1 

49.8 

11.6 

19.5 

0.147 

0.349 

0.660 

2.39 

60 

18.3 

50.8 

11.4 

19.5 

0.145 

0.348 

0.651 

2.39 

90 

17.4 

52.4 

11.7 

18.5 

0.144 

0.345 

0.640 

2.38 

100 

16.4 

51.6 

13.4 

18.6 

0.145 

0.343 

0.629 

2.37 

110 

14.3 

52.2 

14.6 

18.8 

0.144 

0.335 

0.600 

2.34 

120 

13.9 

51.1 

16.6 

18.5 

0.144 

0.327 

0.583 

2.34 
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The  results  show  that  Ps  can  be  polarized  with  an  electric  field  of  a  magnitude  less  than 
120kV/cm[16]. 

Besides  lifetimes,  the  intensity  of  the  Ps  can  be  influenced  by  electric  fields.  In 
an  experiment  run  by  Yasuo  Ito  at  the  Research  Center  for  Nuclear  Science  and 
Technology  at  the  University  of  Tokyo,  Na-22  was  deposited  inside  an  envelope  of  nickel 
foils,  which  was  then  sandwiched  between  two  low  density  polyethylene  sample  sheets. 
By  applying  a  bias  across  the  foils,  the  o-Ps  intensity  not  only  decreased  as  the  electric 
field  strength  increased,  but  the  polarity  of  the  field  did  not  significantly  affect  the 
intensity,  as  shown  in  Figure  6  [reproduced  from  reference  5].  In  this  experiment,  the  o- 
Ps  intensity  is  a  measure  of  the  cross  section  for  o-Ps  formation.  As  the  electric  field 
strength  increases,  the  field  imparts  more  energy  to  the  positron  and  when  it  ionizes  an 
atom  or  molecule,  the  positron  gives  more  energy  to  the  electron.  As  the  energy  of  the 
electron  increases,  the  electron  is  more  likely  to  escape  before  the  positron  can  capture  it 
and  form  o-Ps. 
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Figure  6.  Comparison  of  the  effects  of  positive  and  negative  static  and  alternating  electric  fields  on 
Ps  formation  at  room  temperature. 


22 


3.  Equipment 


The  two  aspects  of  this  research  involve  conducting  PsARS  and  PALS 
measurements.  Each  system  has  different  requirements  for  its  equipment  due  to  the 
particular  aspect  of  the  positron  annihilation  that  the  system  is  measuring.  For  the 
PsARS  system,  the  important  factor  is  being  able  to  capture  the  maximum  amount  of 
photons  from  the  source.  In  PALS,  the  system  must  be  able  to  resolve  the  time  different 
between  the  positron’s  birth  and  death  photons.  These  different  necessities  lead  to  having 
separate  detectors  for  each  system. 

3.1  3y/27  Detector 

3.1.1  Detector  Overview 

To  measure  the  Sy/ly  ratio,  a  positron  annihilation  of  radiation  spectroscopy,  or 
PsARS,  system  is  used.  The  digital  PsARS  system  at  AFIT  consists  of  a  thallium-doped 
Nal  scintillation  ring  that  is  18  cm  long,  has  a  20  cm  diameter,  and  a  thickness  of  3.5  cm. 
Figure  7  shows  the  location  of  the  crystals  within  the  ring,  and  Figure  8  shows  a  picture 
of  the  detector.  The  full  detector  schematic  is  attached  in  Appendix  A.  The  six  Nal(Tl) 
crystals  were  connected  to  six  photomultiplier  tubes,  or  PMTs.  The  PMT’s  are  each 
connected  to  an  ORTEC  model  113  pre-amplifier.  Each  pre-amp  is  connected  to  an  input 
channel  in  an  XIA  DGF4C  revision  E  digital  system,  to  convert  the  current  pulses  into 
digital  data.  Two  cards  are  used  because  each  card  can  only  hold  four  channels.  The 
signals  are  digitized  using  a  12  bit  ADC  at  a  rate  of  500  MHz  [17].  The  XIA  cards  are 
controlled  by  a  Weiner  CC32  CAMAC  controller.  The  signal  from  the  controller  is  sent 
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to  a  computer,  where  it  was  analyzed  using  Igor  Pro,  a  data  analysis,  image  processing. 


and  programming  software  tool. 


Figure  7.  Schematic  of  the  AFIT  ring  detector  made  of  six  Nal(Tl)  scintillation  crystals.  Each  crystal 
is  optically  separated  from  each  other  and  connected  to  its  own  PMT. 


Figure  8.  Picture  of  the  AFIT  Nal(Tl)  ring  detector. 
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The  right  XIA  card  was  set  as  the  master  card  and  the  left  card  was  set  as  the 
terminator  card.  Table  3  shows  the  configuration  of  the  jumpers  on  the  cards  to  achieve 
this  setup. 


Table  3.  Setup  of  the  jumpers  on  the  two  XIA  cards  to  make  the  right  card  the  master  and  left  card 

the  terminator 


Right  card  (s/n:  1230) 

Left  card  (s/n:  1240) 

Jumper  1 

Attached 

Detached 

Jumper  2 

Attached 

Detached 

Jumper  3 

Detached 

Attached 

Jumper  4 

Detached 

Attached 

Jumper  5 

Board  clock 

External 

In  addition  to  the  five  jumpers  listed  above,  all  the  jumpers  to  control  the  input 
impedances  were  set  to  50  Ohms.  To  allow  the  two  cards  to  communicate  to  each  other, 
the  backplane  connector  was  connected  to  the  row  A  of  the  16  pin  connector  on  each 
card,  as  shown  in  Figure  9. 


Backplane  connector 


Figure  9.  Diagram  placement  of  the  backplane  connector  on  the  16  pin  connector  in  a  two  module 
system. 


Furthermore,  a  logic  unit  was  added  to  the  system  so  that  when  one  module  triggered,  it 
forced  the  other  module  to  trigger  as  well.  This  was  accomplished  by  connecting  the 
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busy-out  signal  of  the  DGF4C  cards  to  the  input  of  a  Phillips  Scientific  model  757  mixed 
logic  fan-in/fan-out  unit  and  the  output  of  the  logic  unit  to  the  sync-in  port  on  the  DGF4C 
cards.  To  force  both  modules  to  trigger  on  any  event,  the  coincidence  pattern  was  set  to 
include  all  possible  coincidence  configurations.  This  was  accomplished  using  the 
settings  tab  of  the  DGF  viewer  (all  boxes  were  selected  including  the  “0000”  pattern). 
The  coincidence  window  was  set  to  50  ns,  also  located  in  the  coincidence  pattern 
function.  With  these  settings,  when  any  one  of  the  six  PMT’s  registers  a  hit,  it  sends  out 
a  signal  to  the  logic  unit,  and  the  logic  unit  forces  all  six  PMT’s  to  collect  a  measurement 
for  50  ns. 

A  spectrum  is  collected  from  the  run  tab  of  the  DGF  viewer.  The  user  can  control 
the  number  of  spills,  polling  time,  and  run  time.  A  spill  is  when  the  DGF4C  cards 
transfer  their  data  to  the  computer.  The  cards  can  hold  256  kB  of  list  mode  data  before  its 
buffer  is  full  and  its  memory  is  transferred  to  the  computer  [21].  The  polling  time  is  the 
time  between  each  measurement.  The  run  time  is  how  long  the  data  is  collected  for.  The 
number  of  spills  was  set  high  enough  to  allow  the  run  time  to  decrement  to  zero  before 
the  spills  reached  zero.  The  polling  time  was  set  to  0.05  psec.  At  the  top  of  the  tab  is  a 
drop  down  menu  with  different  types  of  collection  methods.  The  run  type  chosen  was  the 
“energy  and  time  only”  method.  This  method  was  chosen  because  its  output  file  gives 
the  channel  number,  energy,  and  trigger  time  for  every  event  detected.  Furthermore,  the 
energy  and  time  only  aspect  eliminates  excess  data  in  the  output  that  would  increase  the 
file  size.  The  output  data  is  organized  into  columns.  The  first  column  contains  the  event 
number,  the  second  column  contains  the  channel  number  of  the  XIA  signal,  the  third 
column  contains  the  energy  of  the  event  (represented  as  channel  number  of  the  spectra). 
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and  the  fourth  column  identifies  the  time  of  the  event.  All  of  the  events  for  module  one 


are  located  in  the  top  half  of  the  file,  while  the  events  for  module  two  are  located  in  the 


bottom  half  of  the  file.  Figure  10  shows  an  excerpt  from  an  example  output  file. 


Module:  1 

Run  Task:  259 

Run  Start  Time:  0.001634 


Event  No  Channel  No  Energy  Trigger  Time  XIA_PSA  User_PSA 


1 

0 

0 

16336 

1 

1 

0 

16336 

1 

2 

0 

16336 

1 

3 

0 

16336 

2 

0 

0 

36911 

2 

1 

0 

36911 

2 

2 

0 

36911 

2 

3 

0 

36911 

3 

0 

0 

3275 

3 

1 

0 

3275 

3 

2 

0 

3275 

3 

3 

0 

3275 

4 

0 

0 

11104 

4 

1 

0 

11104 

4 

2 

0 

11104 

4 

3 

0 

11104 

5 

0 

0 

19333 

5 

1 

22382 

19333 

5 

2 

3495 

19333 

5 

3 

0 

19333 

6 

0 

0 

51178 

6 

1 

0 

51178 

6 

2 

3778 

51178 

6 

3 

0 

51178 

7 

0 

0 

47619 

7 

1 

0 

47619 

7 

2 

0 

47619 

7 

3 

0 

47619 

8 

0 

0 

27938 

3 

1 

0 

27938 

8 

2 

28466 

27938 

8 

3 

0 

27938 

9 

0 

0 

22325 

9 

1 

38315 

22325 

9 

2 

35158 

22325 

9 

3 

0 

22325 

10 

0 

0 

33794 

Figure  10.  Example  of  an  output  file  from  the  Igor  Pro  system  using  the  energy  and  time  only  run 
type. 


As  can  be  seen  from  the  figure,  there  is  header  data  at  the  top  of  the  file  that  identifies  the 
module  and  run  start  time  of  the  experiment.  The  run  start  time  does  not  report  what  time 
the  experiment  started,  but  whether  the  two  modules  are  synchronized  at  the  start  of  the 
experiment.  If  the  two  modules  are  synced  at  the  start  of  the  data  collection,  the  two  run 
start  times  will  read  0.001634.  This  same  information  exists  in  the  middle  of  the  file  at 
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the  start  of  module  two’s  data.  For  the  Sy/ly  code  to  work  properly,  this  information 
must  be  removed  from  the  output  file. 

3.1.2  Detector  Calibration,  Analysis  Code,  and  Accidental  Coincidence 
Correction 

The  first  step  in  calibrating  the  equipment  involved  determining  the  optimal 
operating  voltage  for  each  PMT.  The  electronics  used  were  an  ORTEC  model  142  pre¬ 
amplifier,  ORTEC  model  556  high  voltage  (HV)  power  supply,  ORTEC  model  671  linear 
amplifier,  ORTEC  model  551  single  channel  analyzer  (SCA),  and  an  ORTEC  model  996 
counter.  Eigure  1 1  shows  the  experimental  setup. 


Pre-ainp  Linear  amp  counter 


Figure  11.  Diagram  of  experimental  setup  to  determine  the  best  operating  voltage  of  each  PMT 

A  Sr-85  source  was  measured  because  its  514  keV  gamma  ray  is  close  to  that  of 
the  511  keV  Ps  annihilation  photo  and  its  activity  was  higher.  To  ensure  the  SCA 
triggered  on  the  514  keV  gamma  ray,  a  Tektronix  TDS5104B  digital  phosphor 
oscilloscope  was  used  to  set  the  correct  cut-off  voltage.  To  minimize  the  background 
radiation,  the  Nal(Tl)  detector  and  subsequent  experiment  was  placed  inside  a  lead  cave, 
reference  Eigure  12.  Eigure  13  and  14  show  background  spectra  collected  outside  and 
inside  the  lead  cave,  using  a  BaEi  detector.  Each  measurement  lasted  60  minutes. 
Comparing  the  two  spectra,  the  lead  cave  suppressed  the  low  energy  background 
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radiation.  Both  figures  show  a  high  energy  peak  at  a  trigger  voltage  of  approximately  0.2 
V.  However,  this  peak  is  a  result  of  how  the  spectra  are  analyzed.  All  of  the  high  energy 
gamma  rays  are  grouped  together  in  this  one  bin,  making  the  appearance  of  high  energy 
background  source  existing.  An  initial  voltage  of  500  V  was  applied  to  the  PMT,  and  the 
counter  recorded  the  number  of  counts  in  10  sec.  The  voltage  was  incremented  until  the 
number  of  counts  reached  a  peak.  The  operating  voltage  of  each  PMT  is  listed  in  Table 
4.  When  looking  at  the  Nal(Tl)  detector  from  the  front,  PMT  1  is  located  at  the  top  and 
the  numbering  system  continues  clockwise,  reference  Figure  15. 


Figure  12.  Picture  of  the  lead  cave  the  vacuum  chamber  system  was  set  up  in.  Next  to  the  lead  cave  is 
the  crate  containing  the  power  suppliers,  pre-amplifiers,  CAMAC  crate,  XIA  cards,  and  vacuum 
sensor  display. 
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Figure  13.  Background  spectrum  collected  using  a  BaF2  detector  located  outside  of  the  lead  cave.  The 
measurement  lasted  60  minutes. 


Figure  14.  Background  spectrum  collected  using  a  BaF2  detector  located  inside  of  the  lead  cave.  The 
measurement  lasted  60  minutes. 


30 


Table  4.  Operating  bias  for  each  Nal(Tl)  PMT 


PMT 

Bias  (V) 

1 

858 

2 

836 

3 

900 

4 

963 

5 

935 

6 

890 

Figure  15.  Numbering  system  used  to  identify  the  PMT’s  when  looking  at  the  front  of  the  Nal(Tl) 
detector. 


Due  to  the  fact  that  only  three  dual  power  supplies  were  available  to  power  the  six 
PMT’s,  the  PMT’s  had  to  be  paired  up  based  on  similar  operating  voltages.  PMT  1  and  2 
were  paired  and  used  831  V,  PMT  3  and  6  were  paired  and  used  890  V,  and  PMT  4  and  5 
were  paired  and  used  930  V.  The  bias  chosen  for  each  pair  could  not  be  an  average  and 
had  to  be  the  lower  of  the  two  optimal  biases.  Once  the  bias  on  the  PMT  increased  past 
its  optimal  value,  the  count  rate  decreased.  The  rest  of  the  calibration  settings  were 
determined  using  the  Igor  Pro  GUI. 

The  filter  and  pulse  shape  analysis  settings  were  chosen  in  order  to  minimize  the 
full  width  at  half  maximum,  FWHM,  and  thus  minimize  the  detector  resolution.  These 
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parameters  were  optimized  by  2"‘^  Lt  Robert  Slaughter  the  previous  year  in  his  thesis 
experiment.  Table  5  and  6  contain  the  values  taken  from  his  thesis: 


Table  5.  3  gamma/2  gamma  experiment  filter  settings  [2:120] 


Trigger  filter 

Energy  filter 

Rise  time  (ps) 

0.1 

2.8 

Flat  top  (ps) 

0.1 

1.2 

Threshold 

40 

N/A 

Table  6.  3  gamma/2  gamma  experiment  pulse  shape  analysis  settings  [2:120] 


Trace  length  (ps) 

5 

Delay  (ps) 

2 

PSA  start  (ps) 

0.025 

PSA  end  (ps) 

4.975 

The  remaining  parameters  were  determined  using  the  built  in  functions  from  the  DGF 
viewer.  The  gains  for  all  the  channels  were  adjusted  so  that  full  1.274  MeV  gamma  ray 
peak  could  be  seen  on  the  DGF  viewer’s  oscilloscope.  If  the  gain  was  set  too  high,  the 
peak  of  the  1.274  MeV  gamma  ray  would  be  cut  off  and  the  signal  would  saturate  the 
electronics.  The  offset  values  were  also  determined  using  the  oscilloscope  function.  The 
offsets  were  first  set  using  the  “adjust  offset  function”,  but  manually  adjusted  for  better 
results.  The  decay  times  of  the  signals  were  set  manually  using  the  “manual  Tau  finder” 
DGF  viewer  function.  The  decay  times  had  to  be  set  individually  to  account  for  the  slight 
differences  between  each  PMT.  Table  7  shows  the  gain,  offset,  and  Tau  settings  for  each 
of  the  channels.  Lastly,  on  the  run  tab  of  the  interface,  the  “simultaneously  start/stop 
modules”,  “synchronize  clocks”,  and  “in  every  new  run”  boxes  were  checked  to  ensure 
the  modules  were  synchronized  and  the  polling  time  was  set  to  0.05  sec. 
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Table  7.  Gain,  offset,  and  Tau  settings  for  the  Nal(Tl)  detector 


Gain  (V/V) 

Offset  (V) 

Tau  (ps) 

Mod  1  eh  1 

0.1640 

1.0000 

68.1501 

Mod  1  eh  2 

0.1640 

1.0000 

64.7665 

Mod  2  eh  0 

0.1640 

1.0348 

67.4105 

Mod  2  eh  1 

0.1640 

1.0590 

63.8633 

Mod  2  eh  2 

0.1640 

1.0552 

74.6553 

Mod  2  eh  3 

0.1640 

1.0301 

66.7308 

The  next  step  of  calibration  involved  collecting  spectra  from  different  isotopes  to 
create  a  calibration  curve  relating  the  channel  number  to  the  energy.  Four  radionuclides 
were  measured:  Co-57,  Cs-137,  Na-22,  and  Sr-85.  Their  respective  gamma  ray  energies 
are:  122  keV,  662  keV,  511  keV,  and  514  keV.  These  sources  were  chosen  to  ensure  the 
calibration  equation  had  a  low,  medium,  and  high  energy  data  point.  Each  source  was 
placed  on  a  holder  developed,  by  2"‘^  Lt  Robert  Slaughter,  to  ensure  the  source  was 
centered  inside  the  detector,  reference  Figure  16. 
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Each  source  was  measured  for  15  minutes  to  collect  an  energy  spectra.  Figure  17 
shows  the  Na-22  energy  spectrum  for  channel  1  of  module  1,  segment  number  1  from  the 
Nal(Tl)  detector.  The  Igor  pro  program  allows  the  user  to  be  able  to  zoom  in  on  the  full 
energy  peak,  FEP,  to  pick  out  its  location.  These  spectra  were  used  to  create  the  six 
energy  calibration  equations. 
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Figure  17.  Uncalibrated  energy  spectrum  for  channel  number  1,  within  module  1,  from  measuring  a 
Na-22  source  for  fifteen  minutes. 

Given  the  location  of  each  FEP  for  the  four  isotopes  and  all  six  channels,  six  calibration 
equations  were  determined.  The  calibration  equations  were  determined  using  a  best  fit 
second  order  polynomial  equation.  The  calibration  coefficients  and  R-squared  value  for 
each  equation  are  located  in  Table  8.  R-squared  is  a  measure  of  how  close  the  best  fit 
equation  fits  the  data.  The  closer  the  value  is  to  1,  the  better  the  fit  is.  The  coefficients 
take  follow  the  format  of: 

Ax^  +Bx  +  C  [15:442]  (19) 


#0 

#1 

#2 

#3 

IH|||||||  1  M 

0  10  20  30  40  50  60 

k 


34 


Table  8.  Coefficients  for  each  Nal(Tl)  detector’s  calibration  equation 


A 

B 

C 

R^ 

Mod  1  eh  1 

-1.27325E-6 

6.95083E-2 

-6.57573E1 

9.99984E-1 

Mod  1  eh  2 

6.82893E-7 

1.95708E-2 

1.62295E1 

9.99994E-1 

Mod  2  eh  0 

-9.42006E-7 

6.32250E-2 

-6.03483E1 

9.99973E-1 

Mod  2  eh  1 

-1.07175E-6 

6.32157E-2 

-6.79067E1 

9.99996E-1 

Mod  2  eh  2 

1.63403E-7 

1.89323E-2 

2.38484E1 

9.99974E-1 

Mod  2  eh  3 

-6.91681E-7 

4.99622E-2 

-6.32222E1 

9.99999E-1 

These  spectra  can  also  be  used  to  determine  the  resolution  of  the  detector.  The  resolution 
is  a  measure  of  how  close  in  energy  two  gamma  rays  can  be  and  also  be  resolved. 
Resolution  is  reported  as  a  percentage,  and  smaller  values  are  better  than  larger  values. 
Resolution,  R,  is  calculated: 

FWHM 

R  =  (20) 

fio 

where  Ho  is  the  location  of  the  FEP  and  FWHM  is  the  full  width  at  half  maximum  of  the 
FEP  [15:115].  The  resolution  of  the  detector  at  122  keV  and  511  keV  for  each  channel 
are  located  in  Table  9. 


Tabie  9.  Resotution  of  the  Nal(Tt)  detector  crystats,  for  each  channet,  at  122  keV  and  511  keV 


122  keV  (%) 

511  keV  (%) 

Mod  1  eh  1 

19.2 

6.9 

Mod  1  eh  2 

12.9 

8.8 

Mod  2  eh  0 

16.9 

7.2 

Mod  2  eh  1 

17.9 

5.7 

Mod  2  eh  2 

16.0 

7.7 

Mod  2  eh  3 

15.8 

5.4 

These  resolution  values  are  important  factors  in  how  the  code  used  to  analyze  the  Igor 
Pro  output  was  developed. 
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A  new  3y/2y  code  was  developed,  using  MATLAB^m,  to  analyze  the  data  because 
the  previous  experimental  setup  where  four  detector  segments  gathered  three  photon 
events  and  two  segments  gathered  two  photon  events  was  inefficient.  The  full  code  can 
be  found  in  Appendix  B.  The  new  PsARS  code  was  developed  by  combining  the  best 
features  of  the  non-coincidence  and  coincidence  methods.  The  code  applied  both 
geometry  and  energy  filters  to  the  data,  and  all  six  segments  of  the  detector  were  used  to 
gather  two  and  three  photon  data.  This  reduced  the  runtime  of  the  experiments  and 
increased  the  detectors  efficiency.  In  the  previous  setup,  PMT  1  and  4  only  triggered  off 
of  2y  annihilation  events  and  the  rest  of  the  PMT’s  only  triggered  off  of  3y  events. 
Although  plenty  of  photons  were  collected,  the  system  was  inefficient  at  collecting 
photons  and  thus  the  measurement  time  was  long.  With  the  new  setup,  all  six  PMT’s 
collected  general  spectra,  and  in  post  processing  the  number  of  and  3y  counts  is 
determined. 

The  new  3yl2y  analysis  code  developed  in  this  research  was  based  on  two  filters, 
geometry  and  energy  filters.  There  is  no  need  for  a  time  filter  in  the  code  because  the 
settings  in  the  code  forced  the  modules  to  synchronize  events  and  take  measurements  at 
the  same  time.  The  code  starts  by  loading  the  data  file  and  initializing  the  variables  that 
record  the  number  of  2y  and  counts.  There  is  only  one  2y  geometry  filter,  which 
counts  the  number  of  events  where  two  opposing  PMT’s  triggered,  but  there  are  three 
different  types  of  geometry  filters.  The  first,  referenced  as  the  120  filter,  keeps  track 
of  3y  counts  that  register  events  on  every  other  PMT,  as  shown  in  Figure  18.  The  second 
fdter,  referenced  as  the  adjacent  fdter,  counts  events  where  two  opposing  PMT’s 
triggered  and  one  adjacent  PMT  that  triggered,  too;  an  example  of  one  such  case  is  shown 
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in  Figure  19.  The  third  filter  is  referenced  as  the  two  count  Sy  filter.  It  looks  for  3y 
photons  that  triggered  non-adjacent  PMT’s.  Figure  20  shows  one  example  of  the 
configuration.  The  case  of  only  two  adjacent  PMT’s  triggering  is  not  accounted  for  in 
the  code  because  this  case  is  most  likely  due  to  Compton  scattering. 


Figure  18.  Diagram  of  the  120  degree  3y  count  filter.  This  variable  records  the  number  of  counts 
that  occur  on  every  other  PMT.  The  black  areas  represent  an  area  that  detected  an  event. 


Figure  19.  Diagram  of  the  adjacent  3y  count  filter.  This  variable  records  the  number  of  counts  that 
occur  on  every  other  PMT.  The  black  areas  represent  an  area  that  detected  an  event. 
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Figure  20.  Diagram  of  the  two  hit  3y  count  filter.  This  variable  records  the  number  of  counts  that 
occur  on  non-adjacent  PMT’s.  The  black  areas  represent  an  area  that  detected  an  event. 


The  next  part  of  the  code  searches  for  the  module  with  the  smaller  number  of 
recorded  events.  Although  the  Igor  Pro  program  synchronized  the  modules  at  the  start  of 
the  experiment,  malfunctions  with  the  programs  software  occasionally  cause  one  module 
not  to  record  an  event  that  the  other  module  records.  This  causes  the  two  XIA  modules  to 
have  differing  numbers  of  events.  To  prevent  the  code  from  crashing  due  to  it  looking 
for  a  non-existent  event,  the  smaller  module  event  total  is  chosen  to  run  the  code.  Even 
though  the  events  from  one  module  to  the  other  may  not  match  up,  coincidence  events 
can  be  matched  based  on  the  XIA  trigger  times.  Following  the  search  for  the  module 
with  the  smaller  number  of  events,  the  energy  calibrations  are  added  in. 

Subsequently,  the  filter  code  compares  the  trigger  times  and  event  number  to 
determine  if  the  events  are  coincident.  The  code  first  compares  the  trigger  times  for  a 
particular  event  to  see  if  the  times  match.  If  the  times  do  not  match,  the  filter  code  uses 
the  event  in  module  one  as  the  reference  event  and  searches  for  the  event  in  the  second 
module  with  the  matching  trigger  time.  Next,  a  second  time  check  is  conducted  to  ensure 
all  six  times  match.  This  second  check  is  necessary  because  occasionally,  a  few  of  the 
trigger  times  within  an  event  will  be  different  from  the  others.  If  this  happens,  or  no 


matching  time  can  be  found,  the  code  moves  on  to  the  next  event.  An  example  of  this  is 


shown  in  Figure  21. 
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Figure  21.  An  example  of  where  times  within  an  event  will  not  match  up  for  the  PsARS  system.  The 
trigger  time  error  occurs  in  event  7.  Channels  0,  1,  and  2  all  have  the  same  trigger  time,  while 
channel  3  has  a  different  time. 


The  rest  of  the  code  looks  for  2y  and  Sy  coincidence  events.  For  each  event  that 
has  had  its  trigger  times  matched  between  modules,  the  code  counts  how  many  channels 
triggered.  Events  with  only  two  or  three  total  hits  were  analyzed.  If  two  hits  were 
counted,  the  event  is  first  tested  by  the  2y  filter.  The  first  filter  is  a  simultaneous 
geometry  and  energy  filter.  It  ehecks  to  see  that  the  triggered  PMT’s  were  opposing,  and 
that  each  photon  was  greater  than  400  keV  and  less  than  600  keV.  The  energy  level 
around  511  keV  was  expanded  to  account  for  the  energy  resolution  of  the  Nal(Tl) 
detector  and  the  energy  calibration.  The  second  filter  checked  to  make  sure  the  sum  of 
the  energies  was  greater  than  one  FEP  and  one  Compton  scatter,  but  less  than  two  FEP’s, 
1022  keV.  The  Compton  scatter  energy  was  set  to  150  keV.  Adding  this  to  the  511  keV 
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FEP,  the  lower  energy  filter  level  was  661  keV.  However,  due  to  the  resolution  of  the 
Nal(Tl)  detector  at  150  keV  (20%)  and  the  lower  limit  of  the  of  the  511  keV  photon  (400 
keV),  the  lower  energy  level  of  filter  was  set  to  580  keV.  The  upper  level  was  set  to 
1124  keV  based  on  the  10%  resolution  at  511  keV.  If  these  energy  and  geometry 
conditions  were  met,  the  event  was  counted  as  a  event. 

If  the  two  hit  event  did  not  pass  the  2y  filter,  it  was  then  tested  by  the  two  count 
Sy  filter.  The  two  count  3y  filter  first  applied  a  different  simultaneous  energy  and 
geometry  filter.  The  filter  checked  to  see  if  the  energies  of  the  two  non-adjacent  hits 
were  greater  than  the  background  radiation,  50  keV,  and  less  than  most  probably  3y 
energy,  380  keV.  Like  before,  the  380  keV  energy  level  was  raised  to  425  keV  to 
account  for  the  detector  resolution.  The  code  then  summed  up  the  energies  of  the 
photons  and  tested  to  see  whether  the  total  was  greater  than  the  sum  of  a  3y  Compton 
scatter  and  one  low  energy  3y  photon  (200  keV),  but  less  than  two  FEP’s  (760  keV).  The 
lower  and  upper  energy  levels  were  expanded  to  allow  for  the  detector  resolution  and  set 
to  160  keV  and  850  keV,  respectively.  If  these  two  conditions  were  met,  the  event  was 
counted  as  a  two  hit  3y  event. 

The  other  two  3y  filters  operate  in  the  same  way  as  the  previous  two.  They  apply 
a  simultaneous  energy  and  geometry  filter,  followed  by  a  second  total  energy  filter.  The 
first  120  degree  filter  checks  that  the  hits  on  every  other  detector  are  greater  than  the 
background,  50  keV,  and  less  than  the  energy  of  the  120  degree  3y  photon,  340  keV.  The 
second  energy  filter,  within  the  120  degree  filter,  then  checks  to  see  if  the  sum  of  the 
energies  is  greater  than  the  sum  of  two  FEP’s  energy  plus  one  low  energy  3y  photon  (730 
keV),  but  less  than  the  sum  of  three  FEP’s  energy  (1022  keV).  With  the  detector 
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resolution  taken  into  account,  the  energy  lower  energy  level  decreased  to  640  keV  and 
the  upper  energy  level  increased  to  1124  keV.  The  second  3y  filter  is  called  the  adjacent 
filter.  The  first  adjacent  filter  tests  has  the  same  50  keV  and  380  keV  energy  levels 
criteria  as  the  two  count  3y  filter,  except  it  checks  that  two  opposing  PMT’s  and  one  of 
the  four  adjacent  PMT’s  triggered.  The  second  filter  within  the  adjacent  filter  checks  to 
see  if  the  sum  of  the  energies  is  greater  than  the  sum  of  two  Compton  scatters  and  one 
low  energy  3y  photon  (350  keV),  but  less  than  three  FEP’s  (1022  keV).  Again,  each 
energy  level  was  adjusted  to  account  for  detector  resolution.  The  resultant  energies  for 
the  lower  and  upper  sum  energy  filter  were,  respectively,  308  keV  and  1124  keV.  Lastly, 
the  code  sums  up  the  total  3y  counts,  and  outputs  the  number  of  2y  and  3y  counts. 

To  correct  for  background  radiation  and  accidental  coincidence,  an  equation  was 
developed  that  related  the  input  count  rate  and  measurement  time.  A  Sr-85  source  was 
measured  and  analyzed  to  develop  the  equation.  The  Sr-58  source  was  chosen  because 
its  gamma-ray  energy  is  close  to  that  of  the  511  annihilation  photon.  Furthermore,  the 
514  keV  gamma-rays  will  not  undergo  Doppler  broadening.  The  511  keV  photon  from 
Na-22  experiences  Doppler  broadening  when  the  annihilating  positron  and  electron  have 
a  momentum  component  along  the  direction  the  photon  was  measured  [22].  The 
accidental  coincidence  counts  arise  from  background  radiation  and  cosmic  rays  triggering 
the  PMT’s  within  the  coincidence  window.  It  addition,  accidental  coincidence  counts 
arise  from  when  a  2y  or  3y  photon  is  coincident  with  background  radiation  or  cosmic 
rays.  When  measuring  the  Sr-85,  no  true  coincidence  events  should  be  happen,  but 
accidental  coincidence  events  will  occur  due  to  the  aforementioned  reasons.  Therefore,  a 
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correction  for  accidental  coincidence  can  be  developed  that  will  be  subtracted  from  the 
raw  3y  and  counts  to  get  the  true  number  of  counts. 

The  correction  was  produced  by  measuring  the  Sr-85  for  one  hour.  The  data  was 
analyzed  using  the  PsARS  code  to  give  the  total  number  of  Sy  and  2y  counts.  The 
formula  relating  the  measurement  time  and  input  count  rate  to  number  of  accidental 
coincidence  events  takes  on  the  form: 

counts  =  xtc  (21) 

where  t  is  the  measurement  time,  c  is  the  input  count  rate,  and  x  is  a  constant  based  on  the 
number  of  accidental  counts  recorded  from  measuring  the  Sr-85  [20].  The  reason  the 
input  count  rate  was  chosen  over  the  activity  of  the  source  was  because  the  Nal(Tl) 
detector  is  very  sensitive  to  even  low  activity  sources  and  will  produce  pulse  pile-up.  For 
example,  when  the  second  capillary  tube  source  was  measured,  the  event  rate  was  2,964.4 
counts/sec,  however,  the  input  rate  into  the  Nal(Tl)  detector  was  964.5  counts/sec. 
Consequently  the  input  count  rate  was  the  true  measurement  of  the  input  into  the 
detector. 

3.2  Lifetime  Measurements 
3.2.1  Igor  Pro  PALS  system 

The  lifetime  measurements  were  conducted  using  PALS  to  measure  the  lifetimes 
of  Ps  formed  at  the  surface  of  the  gold  coated  capillary  tube.  The  method  uses  two  BaF2 
scintillation  detectors  connected  to  a  Tektronics  oscilloscope  to  trigger  coincident  events 
that  occur  within  100  nsec.  Both  of  the  detectors  act  as  the  start  and  the  stop  modules. 
As  soon  as  one  of  the  detectors  is  triggered  by  the  1.274  MeV  gamma-rays,  the  PALS 
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program  tracks  the  amount  of  time  that  passes  until  the  other  detector  registers  one  of  the 
511  keV  annihilation  photons.  If  no  511  keV  photon  is  detected  within  100  ns,  the 
modules  reset  and  wait  for  another  1.274  MeV  gamma-ray  event. 

The  PALS  software  programs  were  developed  at  Eglin  Air  Force  Base,  by  Dr. 
Michael  Lindsay.  The  program  consists  of  four  parts,  the  Digital  PALS  program  file,  the 
PALS  tool  file,  the  VISA  file,  and  the  PALSfit  program.  The  Digital  PALS  program  file 
is  responsible  for  gathering  the  lifetime  data.  At  the  beginning  of  the  code,  there  is  a 
section  that  lists  all  of  the  user  defined  detector  dependent  parameters.  The  AFIT  BaF2 
parameters,  as  shown  in  Table  10,  were  taken  from  the  previous  version  of  the  code, 
DigitalPALS_1.14.ipf,  and  added  to  the  new  code  sent  in  by  Dr.  Lindsay, 
DigitalPALS_1.16.ipf.  The  VISA  file  was  written  to  allow  programs  to  control  virtual 
instrument  software  architecture,  VISA,  compliant  instruments  connected  to  a  computer 
[13].  For  this  experiment,  the  VISA  file  allows  the  Igor  Pro  PALS  program  to  control  and 
communicate  with  the  Tektronix  TDS5104B  digital  phosphor  oscilloscope  function.  The 
tool  file  contains  all  of  the  non- VISA  related  functions  need  for  importing,  exporting,  and 
analyzing  the  raw  PALS  data.  In  addition  to  a  new  version  of  the  digital  PALS  codes,  the 
PALS  tools  code,  DigitalPALSTools_1.6.ipf,  was  updated.  The  output  of  the  tool  file  is 
a  lifetime  spectrum  of  the  PALS  data.  Lastly,  the  spectrum  is  analyzed  using  PALSfit. 
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Table  10.  Igor  Pro  PALS  parameters  for  the  BaF2  detectors 


Ref  voltage  (V) 

-1.5 

Cutoff  Voltage  (V) 

0.4 

Trigger  amp 

0.08 

DP_G  width 

2.6x10'^' 

DP_G  width  range 

0.6x10'^ 

DP_B  baseline 

-0.0025 

DP_V  scale 

0.1 

DP_V  position 

3.8 

dp  t  fit  before 

5.0x10-'' 

dp  t  fit  after 

1.0x10"^ 

Points  per  frame 

500 

Num  frames 

4395 

DP_Ticker 

58.823 

DP_Purl 

0.0021 

DP_Pur2 

-0.0381 

DP_Pur3 

8.7x10'" 

PALSfit  is  a  program  that  uses  a  least  squares  best  fit  routine  to  find  the  Ps 
lifetimes.  PALSfit  operates  in  two  different  modes  to  extract  different  data  from  the 
spectra  [30].  The  first,  called  PositronFit,  extracts  the  lifetimes  and  intensities  from  the 
spectra.  The  second,  called  ResolutionFit,  calculates  the  time  resolution  function  from 
the  spectra  [30].  Both  modes  fit  a  sum  of  decay  exponentials,  convoluted  with  a 
resolution  function,  plus  a  constant  background,  to  the  spectra.  Along  with  the  spectra  to 
be  analyzed,  PALSfit  needs  a  set  of  input  parameters  to  run.  These  parameter  are  best 
guesses  of  the  parameters  to  be  fitted,  and  may  also  be  constraints  to  the  parameters  [30]. 
Once  ran,  the  program  outputs  the  final  estimate  for  all  the  fitting  parameters  and  their 
standard  deviations.  In  addition,  the  program  outputs  three  numbers  that  describe  how 
well  the  parameters  model  the  spectra:  variance  of  the  fit,  chi-square,  and  significance  of 
the  imperfect  model.  The  program  can  also  plot  the  data  and  results  together  [30]. 
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The  first  step  in  calibrating  the  two  BaF2  detectors  involved  determining  the 
optimal  operating  voltage  for  each  PMT.  The  electronics  used  were  an  ORTEC  model 
142  pre- amplifier,  ORTEC  model  556  high  voltage  (HV)  power  supply,  ORTEC  model 
671  linear  amplifier,  ORTEC  model  551  single  channel  analyzer  (SC A),  and  an  ORTEC 
model  996  counter.  The  electronics  were  connected  in  the  same  way  as  for  the 
calibration  of  the  Nal(Tl)  detector,  reference  Eigure  10.  As  before,  a  Sr-85  source  was 
measured.  To  guarantee  the  SCA  triggered  on  the  514  keV  gamma  ray,  a  Tektronix 
TDS5104B  digital  phosphor  oscilloscope  was  employed  to  set  the  correct  cut-off  voltage. 
The  initial  voltages  on  the  PMT’s  were  set  to  -2200  V,  and  the  counter  recorded  the 
number  of  counts  for  10  sec.  The  voltage  was  then  incremented  until  the  number  of 
counts  peaked.  The  operating  biases  for  the  two  detectors  are  located  in  Table  11. 

Table  11.  Operating  biases  of  the  two  BaF2  detectors 


Serial  number 

Bias  (V) 

BA0935 

-2380 

BA0929 

-2460 

After  running  calibrations,  a  setting  on  the  oscilloscope  was  changed,  leading  the 
Igor  Pro  program  to  freeze  as  soon  as  the  Digital  PAES  code  was  executed.  After  months 
of  unsuccessful  troubleshooting,  a  new  system  was  developed  with  the  help  of  Et 
William  Harrell. 

3.2.2  MATLABTM  PALS  system 

The  PAES  system  developed  by  Et  Harrell  was  originally  designed  to  detect  pair 
production  and  Compton  scattered  events  from  high  energy  gamma  rays.  To  be  able  to 
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determine  if  two  events  are  coincident,  Lt  Harrell  developed  a  program,  in  MATLAB^m^ 
to  determine  the  time  between  two  incoming  signals  triggered  in  an  oscilloscope.  The 
setup  involved  the  two  BaFi  detectors  facing  each  other,  with  the  radioactive  source 
located  underneath  the  detectors,  reference  Figure  22.  Unlike  most  PALs  measurements 
where  the  two  detectors  are  each  placed  on  one  side  of  the  radioactive  source,  the  BaF2 
detectors  are  placed  underneath  the  source  because  they  are  more  efficient  in  this 
configuration  .  The  signals  from  the  two  BaFi  detectors  then  lead  into  two  channels  on  a 
Tektronix  DPO7104  digital  phosphor  oscilloscope,  serial  number  B022288.  Through  a 
TekVISA  application,  a  VISA  application  that  allows  the  computer  to  communicate  with 
the  oscilloscope,  the  MATLAB^m  program  is  able  to  communicate  and  control  the 
oscilloscope.  The  instrument  driver  for  the  oscilloscope  was  obtained  from  Captain  Jesse 
Forster,  who  modified  a  driver  that  controls  a  different  model  of  oscilloscope  because  no 
drivers  exist  for  the  DPO7104  model.  A  diagram  of  how  the  electronics  are  connected  is 
shown  in  Figure  23. 


Figure  22.  Setup  of  the  PALS  BaF2  detectors  with  respect  to  the  capillary  tube. 
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Detectors 


Figure  23.  Diagram  of  how  the  electronics  are  connected  together  to  conduct  the  MATLAB™  PALS 
experiment. 

The  first  step  in  creating  the  lifetime  code  involved  taking  the  original  code 
developed  by  Lt  Harrell  and  altering  it  to  measure  radioactive  sources  and  produce 
spectra  in  order  to  calibrate  the  system.  The  initial  block  of  code  sets  up  communication 
with  the  oscilloscope,  and  initializes  all  of  the  oscilloscope  parameters.  Many  of  these 
parameters  were  set  to  match  those  in  the  DigitalPALS_1.16.ipf  PALS  code.  However, 
some  values  could  not  be  met  and  they  were  set  to  be  as  close  as  possible  to  the  Igor  Pro 
PALS  values.  The  oscilloscope  parameters  are  listed  in  Table  12.  The  time  per  division, 
time  per  point,  and  record  length  must  be  set  on  the  oscilloscope  by  hand  before  running 
the  MATLAB^m  program  because  the  program  cannot  set  all  the  needed  parameters. 
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Table  12.  MATLAB^*^  PALS  parameters  for  the  BaF2  detectors 


mV/division 

50 

ns/division 

20 

psec/point 

100 

Record  length 

2000 

Run  type 

High  resolution 

Horizontal  position 

20 

Enhanced  triggering 

off 

Ch  3  trigger  level 

-0.01 

Ch  4  trigger  level 

-0.01 

Once  the  connection  with  the  oscilloscope  is  established,  the  MATLAB^m  code 
starts  collecting  data,  to  calibrate  the  system,  through  the  use  of  the  data_cal_stream.m 
file,  attached  in  Appendix  C.  The  calibration  code  collects  data  through  the  use  of  an 
“AND’  logic  gate.  The  threshold  trigger  level  for  both  channels  is  set  to  -0.01  V  and  the 
condition  for  the  logic  gate  to  trigger  on  was  set  to  false.  This  makes  the  oscilloscope 
trigger  on  any  pulse  with  a  magnitude  above  10  mV.  Next,  the  code  captures  the  pulse’s 
waveform,  and  saves  the  waveform  to  a  file.  Each  time  a  new  pulse  is  detected,  the  new 
waveforms  are  appended  to  the  previous  data  file.  This  data  file  is  then  analyzed  using 
the  cal_scan.m  file,  located  in  Appendix  D.  This  code  loads  the  file,  and  then  loads 
segments  of  data  to  be  analyzed  instead  of  loading  the  entire  data  set  and  then  analyzing 
the  data.  This  greatly  reduces  the  runtime  of  the  code.  The  faster  analysis  time  is 
accomplished  using  the  textscan  function  in  MATLAB^m.  The  textscan  function  loads 
the  data  into  the  cal_func.m  function  file,  located  in  Appendix  E,  file  and  analyzes  the 
segmented  data.  The  cal_func.m  file  outputs  the  results  back  into  the  cal_scan.m  file  and 
the  textscan  function  passes  the  next  group  of  data  to  the  cal_func.m  file. 
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The  cal_func.m  function  file  smoothes  out  the  waveforms  using  a  Savitzky-Golay 
routine  to  fit  the  data  to  a  ninth  order  polynomial  using  thirty-three  waveform  points. 
This  setup  was  chosen  to  smooth  out  the  noise  in  the  data,  but  not  so  much  that  data  is 
lost  from  the  smoothing.  Figure  24  shows  an  example  of  an  unsmoothed  waveform,  and 
Figure  25  shows  the  same  waveform,  but  smoothed.  After  smoothing,  the  code  finds  the 
location  and  pulse  height  of  the  waveform.  Lastly,  the  function  writes  the  pulse  heights 
back  to  the  cal_scan.m  file.  The  scan  file  then  continues  to  load  segment  data  into  the 
function  file  until  all  the  waveforms  have  been  analyzed.  Finally,  the  code  assembles  and 
displays  pulse  height  energy  spectrum  for  each  channel. 


Figure  24.  Unsmoothed  waveform  pulse  from  measuring  a  Na-22  source. 
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Figure  25.  Smoothed  waveform  pulse,  from  measuring  a  Na-22  source,  using  a  Savitzky-Golay  filter. 
The  filter  smoothed  the  data  using  a  ninth  order  polynomial  using  thirty -three  waveform  points. 

From  the  pulse  height  spectra,  calibration  equations  can  be  made  for  the  detectors.  To 
create  the  calibration  equations  for  the  PALS  code,  three  sources  ,Cs-137  (AFIT  source 
number  T-89),  Co-57  (AFIT  source  number  T-85),  and  Sr-85  (AFIT  source  number  T- 
133),  were  measured  for  three  hours  each.  All  three  sources  are  gamma  ray  sources  and 
release,  respectively,  662  keV,  122  keV,  and  514  keV  gamma  rays.  Figure  26  shows  the 
Cs-137  energy  spectrum. 
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Figure  26.  Energy  spectrum  from  measuring  a  Cs-137  source  for  three  hours.  The  662  keV  peak  is 
located  at  0.0482  V. 


From  looking  at  the  Cs-137  spectrum,  it  is  evident  that  there  is  a  bias  in  the  system  that 
causes  certain  trigger  level  bins  to  collect  more  counts  than  others.  The  bias  is  caused 
from  the  rounding  off  of  the  pulse  voltage  in  the  oscilloscope.  Furthermore,  this  problem 
is  aggravated  by  the  system  not  having  enough  resolution  when  the  data  is  digitized.  Due 
to  the  poor  energy  spectra  and  the  large  error  associated  with  identifying  the  FEP  for  each 
isotope,  the  calculated  energy  calibration  equations  did  not  correctly  identify  the  trigger 
voltages  for  the  511  keV  and  1274  keV  photons.  Instead  of  creating  a  calibration 
equation  to  convert  the  trigger  levels  into  energies,  the  PALS  code  ran  directly  off  of  the 
trigger  level  voltages.  The  correct  trigger  levels  were  determined  by  examining  a  Na-22 
source  on  the  oscilloscope  in  infinite  persistence  mode  to  identify  the  voltages.  The  511 
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keV  photon  was  set  to  trigger  on  -35  mV  and  the  1274  keV  photon  on  -95  mV,  but  were 
raised,  respectively,  to  -10  mV  and  -60  mV  to  account  for  the  poor  energy  resolution  of 
the  BaF2  detectors. 

With  the  calibrations  complete,  Lt  Harrell’s  original  code  was  modified  again  to 
create  the  new  PALS  code.  As  before,  the  first  code  sets  the  parameters  on  the 
oscilloscope  and  streams  in  the  data.  This  code  is  called  Stream_data_PALS.m,  and  is 
located  in  Appendix  F.  However,  this  code  does  not  use  the  “AND”  logic  gate  like 
before.  This  code  uses  a  setup  and  hold  trigger  mode.  With  this  setup,  the  first  channel 
acts  as  the  start  module  and  triggers  off  of  pulses  greater  in  magnitude  than  60  mV.  The 
first  trigger  acts  as  the  clock  source  and  once  it  triggers,  off  of  the  1274  keV  photon,  the 
clock  stays  open  either  for  100  ns,  or  until  the  second  channel  triggers.  Off  of  the  511 
keV  photon.  The  second  channel  acts  as  the  stop  module  and  triggers  off  of  pulses 
greater  in  magnitude  than  10  mV.  During  the  time  the  clock  remains  open  the 
oscilloscope  records  the  waveforms  for  both  pulses.  If  a  second  pulse  was  detected 
within  100  ns  that  met  the  energy  criteria,  the  two  waveforms  were  sent  to  the  stream  file 
to  be  saved.  If  a  second  pulse  was  not  detected,  the  oscilloscope  reset  and  waited  for 
another  start  signal  to  begin  recording  data  again. 

Once  the  data  is  collected,  the  data  file  is  analyzed  using  the  PALS_meas.m  file, 
located  in  Appendix  G.  This  file  opens  the  data  file,  and  sends  segmented  data  to  the 
PALS_func.m  file.,  located  in  Appendix  H.  This  file  smoothed  and  applied  energy 
calibrations  to  the  data.  Next,  the  data  was  passed  to  a  second  function  file, 
find_midrise_time.m,  which  is  located  in  Appendix  I.  The  function  found  the  location  in 
the  waveform  where  the  pulse  reaches  50%  of  its  maximum,  then  reported  the  time  for 
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that  location.  The  midrise  time  of  the  pulse,  rather  than  the  time  of  the  pulse  peak,  was 
used  as  the  time  of  the  pulse  because  the  noise  in  the  waveform  changed  the  location  of 
the  pulse  peak.  Next,  the  code  checked  to  see  if  the  difference  between  the  two  midrise 
times  were  less  than  the  100  ns  window.  If  they  were,  an  energy  filter  was  applied  to 
check  to  ensure  the  first  pulse  was  the  1,274  keV  start  pulse  and  the  second  pulse  was  the 
511  keV  stop  pulse.  If  the  energy  criterion  were  met,  the  different  between  the  two  times 
was  reported  as  the  lifetime  of  that  Ps  annihilation  event.  Lastly,  all  the  lifetime  data  was 
passed  back  to  the  PALS_meas.m  file  and  the  data  was  plotted  as  a  lifetime  spectrum. 
This  data  can  then  be  analyzed  using  PALSfit  as  well  to  determine  the  annihilation 
lifetimes  contributions  in  the  PALS  spectrum 

3.3  The  Capillary  Tube 

The  capillary  tubes  were  chosen  for  this  research  because  of  their  ability  to  both 
contain  the  radioactive  Na-22  solution  and  moderate  the  positrons.  The  tubes  were 
custom-made  of  silica  with  an  inner  layer  of  silver.  On  top  of  the  silica  is  a  layer  of 
polyimide.  An  outer  layer  of  gold  was  added  later  by  Lt  Slaughter.  Two  different  tubes 
were  procured,  each  with  a  different  thickness  of  silica.  The  lighter  colored  tube  has  a 
silica  thickness  of  65  pm  and  the  darker  colored  tube  has  a  silica  thickness  of  75  pm, 
reference  Figure  27. 
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Figure  27.  Cut  away  diagram  of  the  custom  made  capillary  tubes  showing  each  material  and  its 
thickness. 


In  Lt  Slaughter’s  experiments,  the  results  were  affected  due  to  the  fact  that  the 
capillary  tube  was  not  centered  inside  the  vacuum  chamber.  The  capillary  tube  was 
attached  on  one  side  to  the  vacuum  chamber  wall,  but  left  unattached  on  the  other  side, 
allowing  the  unattached  end  of  the  capillary  tube  to  bend  downward  under  the  force  of 
gravity.  Figure  28  diagrams  a  device  to  hold  the  free  end  of  the  capillary  tube  up.  The 
ring  was  designed,  by  me,  to  have  a  hollow  center  to  reduce  the  amount  of  mass  escaping 
positrons  could  interact  with.  Fishing  line  was  strung  across  the  vertical  and  horizontal 
holes  to  create  a  square  of  line  in  the  middle  that  will  hold  the  free  end  of  the  capillary 
tube  in  the  center  of  the  vacuum  chamber. 


54 


Figure  28.  Diagram  of  polyethylene  ring  designed  to  hold  the  capillary  tube  in  the  center  of  the 
vacuum  chamber.  Across  the  vertical  and  horizontal  holes  is  strung  fishing  wire  to  hold  the  capillary 
tube  in  place.  The  outer  diameter  is  1.350  inches.  The  inner  diameter  is  0.85  inches.  Each  small  hole 
has  a  diameter  of  1  mm.  Each  of  the  holes  is  located  3  mm  from  edge  of  the  inner  radius,  and  1  mm 
off  the  vertical  and  horizontal  lines  coming  from  the  center.  The  thickness  of  the  ring  is  0.25  inches. 
The  ring  is  constructed  of  polyethylene. 

The  capillary  tube  is  housed  in  a  vacuum  chamber  system  to  reduce  the  amount  of 
atoms  and  molecules  escaping  positrons  could  interact  with.  It  is  attached  at  one  end  to  a 
high  voltage  feed  through  and  the  HV  feed  through  is  connected  to  a  10  inch  ConFlat 
flange  (CF)  nipple,  constructed  by  Kurt  J.  Lesker.  Wires  connecting  the  inside  silver 
layer  and  outside  gold  layer  of  the  capillary  tube  each  connect  to  a  voltage  feed-through 
to  allow  external  biases  to  be  applied  to  the  layers.  After  the  nipple  in  the  system  is  a 
Gamma  Vacuum  lOS  ion  pump,  followed  by  a  valve,  a  Drytel  1025  roughing  pump, 
another  valve,  and  finally  a  Varian  multi-gauge  vacuum,  sensor.  The  roughing  pump  will 
pump  the  pressure  in  the  chamber  down  to  10'^- 10’"^  torr.  Once  the  pressure  reaches  10’4 
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torr,  the  ion  pump  can  be  turned  on  to  reach  pressures  of  10'  torr.  A  tee  can  be 
exchanged  for  a  cross  in  the  middle  of  the  system  to  allow  a  second  experiment  to  be 
going  on  in  the  other  side  of  the  system,  and  the  two  valves  on  both  sides  of  the  roughing 
pump  can  isolate  one  side  of  the  system.  Figure  29  shows  a  picture  of  the  setup. 
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Figure  29.  Experimental  setup  of  the  PsARS  system  with  the  capillary  tube  centered  inside  the 
Nal(Tl)  detector. 


From  Figure  6,  a  majority  of  the  o-Ps  formation  can  be  suppressed  inside  the 
polymer  by  applying  a  50kV/cm  electric  field  across  the  capillary  wall.  Larger  fields  will 
decrease  the  o-Ps  further,  but  with  diminishing  returns.  Above  50  kV/cm,  larger  and 
larger  electric  fields  reduce  the  o-Ps  formation  by  increasingly  smaller  amounts.  To 
determine  the  approximate  biases  needed  to  achieve  a  50  kV/cm,  field,  the  silver  and  gold 
layers  are  modeled  as  a  cylindrical  cathode  and  anode.  Because  the  cathode  and  anode 
are  cylindrical,  and  not  parallel  plates,  the  electric  field  is  not  uniform.  The  electric  field 
from  a  cylindrical  capacitor  calculated: 

E{r)  = - (22) 

r*  Ln(—) 
a 
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where  V  is  the  applied  voltage  between  inner  and  out  cylinders,  a  is  the  inner  cylinder 
radius,  b  is  the  outer  cylinder  radius,  and  r  is  the  location  at  which  the  electric  field  is 
calculated  [15:162].  For  the  capillary  tube,  the  inner  radius  is  3.5x10’"^  m  and  the  outer 
radius  is  4.529x10'"^  m.  Because  the  positrons  are  going  to  be  brought  to  the  outside 
surface  of  the  gold  surface,  r  equals  4.53x10'"^  m.  By  calculating  the  electric  field  at 
varying  voltages  applied  to  the  silver  layer,  the  voltage  needed  to  obtain  an  electric  field 
of  50  kV/cm  will  be  known,  reference  Figure  30. 


Voltage  (V) 

Figure  30.  Calculation  of  the  electric  field  at  the  top  of  the  gold  surface  on  the  capillary  tube  with 
varying  voltages  applied  to  the  silver  layer. 


From  Figure  30,  an  electric  field  of  50  kV/cm  will  be  obtained  by  applying  a  bias  of 
approximately  600  V  to  the  silver  layer.  However,  this  value  is  only  an  estimate  of  the 
true  value  because  the  silver  and  gold  layers  are  not  perfect  cylinders  and  the  gold  layer  is 
not  100  nm  thick  over  the  whole  capillary  tube.  Furthermore,  the  calculation  does  not 
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account  for  the  dielectric  in-between  the  silver  and  gold  layers.  This  calculation  is  only 
used  to  estimate  the  approximate  bias  needed  on  the  silver  layer. 

The  gold  was  added  to  the  capillary  tube,  by  Lt  Slaughter,  using  the  torr  Electron 
Beam  Evaporation  System  (E-Beam)  [2:65].  To  determine  the  smoothness  of  the  gold 
layer  on  the  capillary  tube,  the  tube  was  scanned  using  an  atomic  force  microscope,  or 
AEM.  The  AEM  works  through  the  use  of  a  cantilever  with  one  end  having  a  sharp  tip. 
The  tip  of  the  cantilever  was  made  of  diamond  with  a  diameter  of  100  nm  [2:73-5].  The 
AEM  was  used  to  scan  a  2000  nm  region  on  the  capillary  tube  that  was  previous 
determined  to  be  smooth  with  no  bumps.  In  previous  scans,  different  areas  of  the 
capillary  tube  had  bumps  as  large  as  40  nm  on  the  surface  [2:73-5].  Eigure  31  shows  the 
change  in  surface  roughness. 
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Figure  31.  Plot  of  surface  roughness  over  a  2000  nm  smooth  region  of  the  capillary  tube  gathered 
using  an  AFM  [2:75]. 
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Each  typical  feature,  or  peak,  size  is  about  100  nm  in  the  lateral  dimension  and  10  nm  in 
height.  Over  the  2000  nm  region,  the  average  change  in  surface  thickness  never  exceeds 
20  nm.  Although  the  surface  roughness  makes  estimating  the  electric  field  difficult,  it 
will  prove  useful  once  the  polarizing  electric  field  is  applied  [2:73-75]. 

To  polarize  the  Ps  on  the  surface  of  the  gold,  an  electric  field  of  approximately 
120  kV/cm  will  be  needed  [16].  To  calculate  the  needed  bias,  equation  22  was  used.  For 
this  case,  the  gold  layer  is  the  inner  cylinder  and  the  vacuum  chamber  wall  is  the  outer 
cylinder  because  it  is  the  location  of  the  ground.  R  equals  4.53x10'"^  m,  “a”  equals 
4.53x10'  m,  and  “b”  equals  1.715x10'  m.  By  calculating  the  electric  field  at  varying 
biases.  Figure  32  was  obtained: 
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Figure  32.  Calculation  of  the  electric  field  at  the  top  of  the  gold  surface  on  the  capillary  tube  with 
varying  voltages  applied  to  the  gold  layer. 
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As  can  be  seen  from  the  graph,  the  values  of  the  electric  field,  obtained  from  applying  the 
voltages,  to  a  different  loeation,  are  one  to  two  orders  of  magnitude  lower  than  the 
previous  ealculation.  This  is  beeause  as  the  location  of  the  inner  cylinder  ehanged  from 
the  silver  to  the  gold  layer  and  the  radius  of  the  layer  inereased.  The  eleetrie  field  is 
highly  dependent  on  the  radii  of  the  eylinders  and  the  radius  of  the  outer  layer  inereased 
by  two  orders  of  magnitude.  To  aehieve  an  eleetrie  field  of  120  kV/em,  a  bias  of 
approximately  20  kV  is  needed.  This  eannot  be  achieved  with  the  eurrent  power  supplies 
in  the  laboratory.  However,  beeause  the  gold  eoating  on  the  eapillary  tube  is  not 
perfeetly  smooth,  the  differing  heights  of  the  surfaee  may  ereate  localized  electric  fields 
on  the  surfaee  that  will  reaeh  the  needed  120  kV/em  [29].  These  high  loealized  eleetrie 
fields  will  be  used  to  polarize  the  Ps,  increasing  the  annihilation  lifetimes.  Figures  33 
diagrams  how  the  biases  were  applied  to  the  inner  and  outer  metallie  layers  of  the 
eapillary  tube. 


Figure  33.  Diagram  showing  the  location  of  where  the  inner  and  outer  voltages  are  applied  to  the 
capillary  tube. 
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From  the  previous  year,  to  deposit  the  Na-22  on  the  inner  wall  of  the  tube,  Lt 

22 

Slaughter  used  an  evaporative  method.  Using  a  22  gauge  syringe  radioactive  NaCl  was 
deposited  inside  the  tube.  Due  to  the  long  time  that  was  need  for  the  liquid  to  evaporate, 
the  tube  was  heated  using  a  hot  plate.  0.26  pCi  of  Na-22  was  intended  to  be  deposited 
inside  the  tube,  however,  much  of  the  NaCl  escaped  from  the  capillary.  Lt  Slaughter 
determined  an  activity  of  only  0.10  pCi  remained  inside  the  tube.  Next,  he  used  the  torr 
Electron  Beam  Evaporation  System  (E-Beam)  to  deposit  the  gold  on  the  surface  of  the 
capillary  tube.  Eastly,  he  attached  the  capillary  tube  to  the  inside  of  the  vacuum  tube 
[2:112].  Through  the  use  of  multiple  voltage  feed-through’s,  Lt  Slaughter  demonstrated 
that  electrical  contacts  could  be  made  both  in  the  inside  and  outside  of  the  capillary  tube 
to  produce  a  voltage  difference  across  the  wall.  He  used  a  22  gauge  electrical  wire  to 
make  the  first  connection  between  the  silver  on  the  inside  of  the  tube  and  one  electric 
feed-through.  The  second  connection  was  made  by  wrapping  the  gold  tube  and  a  22 
gauge  wire  with  a  thin  electrical  wire,  and  soldering  the  22  gauge  wire  to  a  second 
voltage  feed-through  [2:112]. 

During  the  deposition  of  the  Na-22  into  the  first  capillary  tube,  more  than  half  of 
the  Na-22  escaped  the  tube.  Lt  Slaughter  made  the  capillary  tube  only  have  contact  with 
the  hot  plate  on  one  side  of  the  tube.  This  caused  a  pressure  build-up  inside  the  tube  and 
ejected  the  Na-22  out.  A  second  capillary  tube  was  created  this  year  to  test  a  new  method 
of  Na-22  deposition.  85  pL,  approximately  0.3  pCi,  of  the  radioactive  Na-22  solution 
was  injected  into  a  new  capillary  tube  using  a  22  gauge  syringe.  The  tube  was  then 
placed  on  two  stands  that  rested  on  a  hot  plate.  Aluminum  foil  was  placed  on  the  hot 
plate  to  prevent  the  hot  plate  from  getting  contaminated.  The  foil  was  also  used  to  create 
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two  contacts  on  which  the  capillary  tube  would  rest.  The  heat  would  travel  from  the  hot 
plate,  through  the  two  metal  contacts,  and  heat  both  ends  of  the  capillary  tube.  By 
heating  both  ends  of  the  tube  instead  of  just  one,  it  is  theorized  that  the  pressure  will  not 
build  up  inside  the  tube  and  the  NaCl  solution  will  not  be  expelled  from  the  tube. 

3.4  Copper  Standard  Source  and  3y/2y  Correction  Factor 

To  correct  for  the  different  efficiencies  at  which  the  Nal(Tl)  detector  captures  a 
photon  versus  a  Sy  photon,  a  photoelectric  absorption  correction  is  needed.  Due  to  its 
lower  energy,  a  3y  photon  is  easier  to  capture  than  a  2y  photon,  but  it  is  harder  to  capture 
all  of  the  three  3y  photons  than  the  two  2y  photons.  To  account  for  the  difference  in 
requiring  the  capture  of  both  2y  photons  but  not  all  three  3y  photons,  a  geometry 
correction  is  also  needed.  In  addition,  the  3y  photon  peak  lies  on  top  of  the  511  keV 
Compton  continuum,  and  because  they  are  not  emitted  collinear,  they  are  harder  to 
identify  than  the  2y  photons.  This  is  corrected  by  adding  in  an  energy  correction  factor. 
However,  all  of  the  correction  factors  can  be  combined  into  one  correction  factor,  the 
3y/2y  correction  factor. 

The  3y/2y  correction  factor  was  calculated  by  measuring  a  standard  source  with  a 
known  3y/2y  ratio.  Single  crystal  copper  is  a  metal  with  a  known  3y/2y  ratio,  1/372.  The 
source  was  created  by  heating  a  copper  plate  on  a  hot  plate.  Slowly,  85  pL,  0.3  pCi,  of 
the  liquid  radioactive  Na-22  was  deposited  onto  the  copper.  The  heat  from  the  hot  plate 
evaporated  the  liquid.  Next,  a  second  copper  plate  was  placed  onto  of  the  first  to 
sandwich  the  Na-22  between  the  copper.  This  was  then  sealed  inside  tape  to  prevent  air 
from  interacting  with  the  copper.  By  collecting  the  spectra  from  the  copper  source  and 
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calculating  the  experimental  3yl2y  ratio,  the  known  ratio  is  divided  by  the  experimental 
ratio  to  determine  the  correction  factor.  Figure  34  shows  what  the  standard  source  looks 
like.  Given  the  accidental  coincidence  correction  and  the  3yl2y  correction  factor,  the 
3y/2y  ratio  can  be  calculated.  The  raw  number  of  3y  and  2y  counts  is  found  from 
analyzing  the  data.  After  the  background  is  subtracted,  the  counts  are  divided  by  the 
2y  counts  and  multiplied  by  the  3yl2y  correction  factor  to  determine  the  actual  3yl2y 
ratio. 


Figure  34.  Picture  of  the  copper  standard  source  used  to  calculate  the  3y/2y  correction  factor. 

This  section  of  the  document  described  the  key  differences  between  the  PsARS 
and  PALS  systems.  The  PsARS  system  used  a  Nal(Tl)  ring  detector  to  maximize  the 
amount  of  photons  collected  while  the  PALS  system  used  fast  BaF2  detectors  to  be  able 
to  resolve  the  time  difference  between  the  1,274  and  511  keV  photons.  In  addition,  a 
standard  source  and  accidental  coincidence  equations  had  to  be  created  to  calibrate  the 
PsARS  system.  The  next  section  will  discuss  the  results  from  measuring  the  copper 
standard  source,  Sr-85  source,  and  capillary  tube. 
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4.  Results  and  Analysis 


Once  each  system  was  setup  and  calibrated,  initial  data  runs  had  to  be  completed 
before  analysis  of  the  capillary  tube  could  begin.  These  include  gathering  information  to 
compute  the  accidental  coincidence  correction  counts  and  the  Sy/ly  correction  factor. 
This  information  was  then  used  to  conduct  the  3y/2y  ratio  and  lifetime  analysis  of  the 
capillary  tube. 


4.1  Accidental  Coincidence  Correction 


The  accidental  coincidence  equation  coefficients  were  calculated  by  measuring 
and  analyzing  a  Sr-85  source.  The  Sr-85  source  was  measured  for  one  hour  and  analyzed 
using  the  3y/2y  code.  By  entering  in  the  Sr-85  input  count  rate,  total  number  of  3y  and  2y 
counts,  and  measurement  time  into  equation  21,  the  coefficients  for  the  3y  and  2y 
equations  were  calculated.  Given  that  the  input  count  rate  from  the  Sr-85  was  3200 
counts/sec,  the  measurement  time  was  3600  sec,  and  the  3y  and  2y  accidental  counts, 
respectfully,  equaled  5  and  35,  the  3y  and  2y  accidental  count  equations  are: 


4.3403x10-^  ct 
=3.0382xl0-^ct 


(23) 


where  c  in  the  input  count  rate  and  t  is  the  measurement  time.  Three  total  measurement 
times  were  used:  2.0  hours,  3.0  hours,  and  3  hours  and  15  minutes,  to  measure  the 
capillary  tube.  With  an  input  count  rate  of  220  counts/sec,  the  accidental  coincidence 
counts  for  the  three  times  are  listed  in  Table  13. 
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Table  13.  3y  and  Ij  accidental  coincidence  counts  for  three  different  total  measurement  times  of  the 
_  capillary  tube. _ _ 


Time  (hrs) 

3y  counts 

2y  counts 

2 

1 

5 

3 

1 

7 

3.25 

1 

8 

Each  of  the  accidental  coincidence  counts  is  low  due  to  the  low  activity  of  the  capillary 
tube.  For  the  3y  counts,  the  2  hour  count  estimate  was  rounded  up  and  the  3  and  3.25 
hour  counts  were  rounded  down.  Differences  in  the  3y  counts  will  only  be  seen  when 
measuring  for  longer  periods  of  time  or  on  higher  activity  sources.  Due  to  the  low  activity 
of  the  source  and  short  measurement  time,  the  error  associated  with  each  count, 
especially  the  3y  counts,  are  large.  When  calculating  the  3j/2j  ratio,  each  of  these 
accidental  coincidence  counts  will  be  subtracted  from  the  raw  counts  to  achieve  the  true 
number  of  counts.  The  low  number  of  accidental  counts  is  a  positive  confirmation  that 
the  criteria  for  determining  if  an  event  was  a  2y  or  3y  event  were  successful.  However, 
the  3y  counts  have  more  error  than  the  2y  counts  because  of  the  low  numbers.  For  two 
hours,  the  calculated  accidental  3y  count  was  less  than  one,  but  rounded  up  to  one.  For 
the  other  two  times,  the  3y  counts  were  calculated  to  be  greater  than  one,  but  were 
subsequently  rounded  down. 

4.2  3y/2y  Correction  Factor 

With  the  ability  to  subtract  out  the  background  radiation,  the  copper  standard 
source  was  measured  using  the  Nal(Tl)  detector  to  create  the  3j/2j  factor.  The  copper 
source  was  measured  for  three  hours  and  fifteen  minutes.  The  data  was  analyzed  using 
the  PsARS  code  to  calculate  the  total  number  of  3y  and  2y  counts,  respectively,  to  be  546 
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and  42395.  Thirteen  separate  measurements  were  conducted  to  obtain  these  sums.  The 


error  in  the  counts  was  calculated  by  dividing  their  standard  deviations  by  the  square  root 
of  the  number  of  measurements.  Dividing  the  number  of  Sy  counts  to  2y  counts  yielded 
an  experimental  3y/2y  ratio  of  0.012858.  Dividing  the  actual  ratio,  1/372,  by  the 
experimental  ratio  yields  the  correction  factor,  0.2090707.  Equation  24  shows  how  the 
correction  factor  is  calculated. 


_  372 
3 

exp- 


(24) 


3 

where  c  is  the  correction  factor  and  exp—  is  the  experimental  3y/2y  ratio.  The  errors  in 
the  number  of  counts  were  then  used  to  calculate  the  error  in  the  correction  factor,  Oc. 


1 


a  3 
exp— 

I*.  2 


exp- 


(25) 


Where  cr  3  is  the  error  in  the  experimental  3y/2y  ratio  and  p  is  the  power  of  the 

exp— 

2 

correction  factor.  Because  the  actual  3y/2y  ratio  is  multiplied  by  the  experimental  ratio 
raised  to  the  power  of  negative  1,  p  equals  the  absolute  value  of  negative  one.  The  error 
in  the  correction  factor  was  calculated  to  be  0.00024,  or  0.115%  of  the  value  of  the 
correction  factor. 

In  addition  to  the  correction  factor,  the  activity  of  the  copper  standard  source  was 
determined  experimentally  as  well.  The  source  was  measured  using  a  high  purity 
germanium  detector,  serial  number  10013369.  The  copper  source  and  a  standard  Na-22 
source  were  placed  on  top  of  the  detector.  By  measuring  the  two  sources  together,  the 
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511  keV  FEP  from  the  energy  spectrum  was  the  result  of  the  photons  from  both  sources. 
Because  the  activity  of  the  Na-22  standard  source  is  known,  the  source  can  be  measured 
using  the  HPGe  detector  separately,  and  the  activity  of  the  copper  source  can  be 
calculated: 

4  ActivityC^Na)  ^  ^  ,  22^,^ 

Activity = -  22.2  4  counts{Cu  +  No)  (26) 

counts{  Na) 

where  the  number  of  Na-22  counts  will  be  determined  from  measuring  the  Na-22  only 
and  the  Cu  plus  Na-22  counts  will  be  calculated  from  measuring  the  two  sources 
together.  After  measuring  the  Na-22  source  for  five  minutes,  1.46x10^  counts  were 
recorded  under  the  511  keV  FEP.  From  measuring  the  two  sources  together,  for  five 
minutes  using  the  same  geometry,  1.48x10^  total  counts  were  recorded.  Subtracting  the 
Na-22  standard  counts  from  the  total  counts  resulted  in  2.0x10"^  copper  standard  source 
counts.  The  activity  of  the  Na-22  was  53.96  pCi  on  1  Jan  05.  After  accounting  for 
decay,  its  current  activity  is  10.824  pCi  (10  Jan  11).  Plugging  these  values  into  equation 
26,  the  activity  of  the  copper  standard  source  is  approximately  0.15  pCi  (10  Jan  1 1).  The 
desired  0.3  pCi  was  not  deposited  due  to  the  error  associated  with  extracting  the  Na-22 
solution  using  a  syringe. 

4.3  3y/27  Ratio  Results 

The  first  3yl2y  ratio  measurements  were  conducted  by  applying  a  voltage  bias  on 
the  inner  silver  layer  with  the  ground  set  to  the  vacuum  chamber  wall.  All  of  the 
measurements  were  conducted  under  a  vacuum  pressure  of  1.36x10'  torr.  The  bias  was 
chosen  to  be  positive  on  the  inner  silver  layer  because  the  positrons  are  positively 
charged.  Using  a  positive  bias,  the  electric  field  will  push  positrons  in  the  direction  from 
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the  silver  layer  to  the  gold  layer.  This  tends  to  force  the  positrons  to  travel  from  the  inner 
portions  of  the  capillary  tube  towards  the  gold  layer.  An  initial  ratio  baseline  was 
conducted  by  applying  no  voltage,  and  then  data  was  collected  with  the  voltage 
incremented  up  to  1  kV,  in  steps  of  200  V.  However,  data  was  collected  at  addition 
biases  to  fill  in  the  gaps  of  the  data.  Each  data  point  was  initially  measured  eight  times  for 
fifteen  minutes,  or  for  a  total  of  two  hours.  However,  to  decrease  the  error  in  the 
measurements  at  the  lower  voltages  to  better  see  what  was  happening,  an  additional  four 
measurements  were  conducted.  These  short  measurement  times  were  chosen  to  minimize 
the  output  data  file  sizes.  For  larger  run  times,  the  output  files  were  too  large  to  open  and 
the  file  information  located  at  the  top  and  middle  of  the  files  could  not  be  deleted.  The 
following  results  were  obtained: 


Figure  35.  3y/2y  ratio  results  by  applying  a  voltage  bias  on  the  inner  silver  layer  of  the  capillary  tube. 
The  ground  was  established  at  the  vacuum  chamber  wall. 
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The  results  show  that  the  Sy/ly  ratio  remains  fairly  constant  up  to  300  V,  then 
rises  at  350  V,  decreases  at  450  V,  and  then  above  450  V  increases  again.  The  change  in 
the  3y/2y  ratio  between  300  and  450  V  show  that  applied  voltages  are  influencing  the 
3y/2y  ratio,  but  the  constant  ratio  before  300  V  show  that  the  ratio  is  changing,  but 
something  is  interfering  with  the  measurements.  However,  as  the  bias  continued  to 
increase  above  450  V,  the  results  show  that  the  bias  was  strong  enough  to  eject  the 
positrons  out  of  the  capillary  tube.  At  high  voltages,  the  electric  field  was  able  to  impart 
enough  energy  to  the  positrons  so  that  they  could  overcome  the  potential  barrier  and  the 
gold  surface,  and  escape.  Once  escaped,  they  either  annihilated  with  molecules  or  atoms 
in  the  air,  or  at  the  vacuum  chamber  wall.  The  rise  in  the  3y/2y  ratio  was  because  the 
gamma  rays  were  being  measured  off  center.  When  the  2y  photons  were  produced  at  the 
center  of  the  detector,  they  were  released  180  degrees  apart  and  detected  by  opposing 
PMT’s.  However,  when  the  positron  annihilated  off  center,  the  2y  photons  were  not 
detected  by  opposing  PMT’s;  they  were  mistaken  for  3y  counts  and  falsely  increased  the 
3y/2y  ratio. 

The  interference  problem  at  lower  voltages  arose  due  to  the  fact  that  the  ground 
was  applied  to  the  vacuum  chamber  wall  and  not  the  outer  gold  layer.  With  a  bias 
applied  to  the  silver  layer,  and  the  gold  layer  left  ungrounded,  the  gold  layer  acted  as  a 
shield  and  prevented  the  voltage  from  dropping  across  the  capillary  tube.  Instead  of  the 
voltage  drop  across  the  capillary  tube,  the  bias  on  the  silver  layer  charged  the  gold  layer, 
and  the  majority  of  the  voltage  drop  occurred  from  the  gold  layer  to  the  vacuum  chamber 
wall.  A  small  experiment  was  conducted  to  test  this  theory.  The  inner  applied  voltage  on 
the  silver  layer  was  set  to  50  V  and  the  gold  layer  was  left  ungrounded.  The  voltage  drop 
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across  the  two  layers  was  .001  V,  but  the  voltage  drop  from  the  silver  layer  to  the  vacuum 
chamber  wall  was  50  V.  When  the  gold  layer  was  established  as  ground,  the  voltage  drop 
across  the  silver  and  gold  layer  increased  to  45  V,  proving  that  with  the  ground 
established  at  the  gold  layer,  the  voltage  drop  across  the  capillary  tube  is  known.  The  45 
V  drop  across  the  silver  and  gold  layer  is  a  best  estimate  because  the  multi-meter  used  to 
measured  the  voltage  does  not  have  an  impedance  high  enough  to  limit  the  current  it 
draws  from  the  circuit,  thus  lowering  the  readout.  To  better  see  what  was  happening  with 
the  Sy/ly  ratio  at  low  applied  voltages,  a  second  experiment  on  the  capillary  tube  was 
conducted  with  the  ground  established  at  the  gold  layer. 

By  establishing  the  ground  at  the  gold  layer,  the  voltage  drop  across  the  capillary 
tube  was  known  and  the  low  voltage  3y/2y  measurements  were  repeated.  Like  the 
previous  experiment,  the  measurements  were  conducted  under  vacuum,  with  a  pressure 
of  1.36x10'  torr.  Each  data  point  was  obtained  by  measuring  the  capillary  tube  for  a 
total  of  three  hours  made  from  twelve  separate  measurements.  The  new  3y/2y  ratios  are: 
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Applied  Bias  on  Silver  Layer(V) 

Figure  36.  3y/2y  ratio  results  by  applyiug  a  voltage  bias  on  the  inner  silver  layer  of  the  capillary  tube. 
The  ground  was  established  at  the  outer  gold  layer. 

The  results  clearly  show  that  the  3yl2y  ratio  rises,  then  falls.  The  peak  3yl2y  ratio  occurs 
at  75  V.  This  peak  is  the  voltage  that  gives  the  positrons  enough  energy  to  lie  at  the  top 
of  the  Ore  gap  as  they  traverse  through  the  capillary  tube.  Beyond  this  voltage,  too  much 
energy  is  given  the  positron  and  as  it  ionizes  an  atom  during  its  travel,  the  positron 
imparts  too  much  energy  to  the  electron  and  the  electron  escapes  before  the  positron  can 
capture  it  and  form  Ps.  To  determine  whether  the  change  in  the  Sy/Sy  ratio  is  influenced 
more  by  the  number  of  3y  or  2y  counts,  the  number  of  counts  are  graphed  in  Figures  37 
and  38.  The  error  data  points  were  calculated  by  taking  the  square-root  of  the  number  of 
counts. 
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Figure  37.  Number  of  2y  couuts  from  the  3y/2y  ratio  experimeut  with  the  grouud  applied  to  the  gold 
layer. 


Figure  38.  Number  of  3y  couuts  from  the  3y/2y  ratio  experimeut  with  the  ground  applied  to  the  gold 
layer. 
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In  Figure  37,  there  is  no  distinct  pattern  in  the  data,  the  data  is  approximately  constant 
when  the  error  bars  are  expanded  from  1  sigma  to  3  sigma.  As  the  3yl2y  ratio  peaks,  then 
decreases,  the  number  of  counts  does  not  peak  and  then  decrease.  In  Figure  38,  the 
data  for  the  counts  follow  the  same  rise  and  fall  pattern  of  the  3yl2y  ratio  data.  This 
implies  that  the  number  of  counts  influences  the  3yl2y  ratio,  and  not  the  number  of  2y 
counts.  The  decrease  in  o-Ps  counts  is  the  result  of  either  the  electric  field  causing  the  o- 
Ps  to  not  form,  or  the  electric  field  is  increasing  the  rate  of  o-Ps  pick  off  annihilation. 
However,  this  determination  cannot  be  made  until  lifetime  analysis  is  conducted  on  the 
data. 

The  voltage  drop  of  75  V  was  used  in  the  subsequent  outer  applied  bias 
experiment.  In  the  second  experiment,  two  biases  were  used,  one  on  the  silver  and  one 
on  the  gold  layer.  The  first  bias  on  the  silver  layer  will  still  push  the  positions  to  the 
surface  of  the  capillary  tube,  and  the  second  bias  will  create  an  electric  field  that  will 
polarize  the  Ps  on  the  gold  surface.  However,  instead  of  applying  a  75  V  bias  to  the 
silver  layer,  the  inner  and  outer  biases  were  adjusted  to  achieve  a  drop  of  75  V  from  the 
silver  to  the  gold  layer.  The  grounds  for  each  bias  were  established  at  the  vacuum 
chamber  wall.  As  before,  the  thirteen  measurements,  fifteen  minutes  each,  were 
conducted  to  gather  each  data  point  and  the  measurements  were  conducted  under  a 
vacuum  of  1.36x10'  torr.  The  3y/2y  ratios  with  an  inner  and  outer  field  applied  to  the 
capillary  tube  are  shown  in  Figure  39. 
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Figure  39.  3y/2y  ratio  results  by  applying  a  voltage  bias  on  tbe  inner  silver  layer  and  outer  gold  layer 
of  tbe  capillary  tube.  Tbe  inner  bias  was  adjusted  so  that  tbe  voltage  drop  from  tbe  silver  to  gold 
layer  was  always  75  V.  Tbe  grounds  were  both  established  at  tbe  vacuum  chamber  wall. 


The  graph  indicates  that  as  the  two  biases  are  increased,  the  Sy/Sy  ratio  decreases,  then 
starts  to  rise  and  remains  fairly  constant.  It  appears  that  the  outer,  polarizing,  electric 
field  is  influencing  the  positrons  before  Ps  is  formed.  The  3y/2y  ratio  decreases  because 
the  outer  electric  field  is  increasing  the  energy  of  the  positron,  pushing  it  above  the  Ore 
gap.  As  the  outer  bias  is  increased,  the  field  imparts  so  much  energy  to  the  positrons  that 
they  are  ejected  from  the  capillary  tube.  As  before,  when  they  annihilate,  because  the 
annihilations  are  off  center,  the  2y  photons  are  mistaken  for  3y  photons,  and  incorrectly 
increase  the  3y/2y  ratio.  In  addition,  the  belief  that  the  positrons  are  being  ejected  from 
the  capillary  tube  is  supported  by  the  fact  that  the  3y/2y  ratio  gradually  increases  as  the 
annihilation  occurs  more  off-center.  If  the  rise  in  the  3y/2y  ratio  was  due  to  a  change  in 
lifetime,  the  ratio  jump  would  occur  suddenly,  and  not  gradually.  This  conclusion  can  be 
confirmed  by  collecting  and  analyzing  lifetime  spectra  from  the  capillary  tube. 
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4.4  PALS  Analysis 


Due  to  the  poor  counting  statistics  achieved  by  Lt  Slaughter,  questions  arose  as  to 
whether  the  digital  PALS  system  was  working  correctly.  To  test  this,  a  Na-22  source  was 
measured  for  1200  sec  using  the  Digital? ALS_1.9.ipf  code,  and  then  analyzed  using  the 
Digital? ALSTools_l. 3. ipf  code.  The  refined  data  was  sent  to  Dr.  Lindsay,  and  he 
determined  that  the  system  was  functioning  correctly.  The  AFIT  PALS  system  was  then 
updated  using  new  code  provided  by  Dr.  Lindsay.  However,  after  the  newer  versions  of 
the  PALS  codes  were  updated  with  the  parameters  for  the  AFIT  BaF2  detectors,  the 
system  no  longer  worked  properly. 

Each  time  the  PALS  code  was  executed,  the  Igor  Pro  program  froze.  All  of  the 
previous  versions  of  the  PALS  code  were  tested  as  well,  but  they  too  froze  up.  Next,  the 
oscilloscope  was  restored  to  its  default  settings,  but  the  PALS  program  continued  to 
freeze.  The  Igor  Pro  program  was  then  uninstalled  and  reinstalled  to  the  oscilloscope,  but 
that  did  not  work.  After  that,  the  entire  PALS  system  was  loaded  onto  a  different 
oscilloscope,  of  the  same  model,  but  the  program  still  did  not  run.  After  months  of 
unsuccessfully  trying  to  correct  the  problem,  no  solution  was  found.  To  be  able  to 
conduct  the  lifetime  analysis,  a  system  used  to  conduct  positron  spectroscopy  was  altered 
to  measure  lifetimes. 

The  MATLAB^m  PALS  program  proved  to  have  problems  as  well.  Two  main 
problems  existed  with  the  MATLAB^m  program:  calibration  and  activity.  When 
measuring  the  standard  Sr-85,  Co-57,  and  Cs-13  sources,  the  energy  spectra  produced 
were  horrible.  As  can  be  seen  in  Figure  25,  identifying  the  662  keV  peak  in  the  Cs-137 
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spectra  is  very  difficult.  In  the  spectra  from  the  other  sources,  it  was  even  harder  to 
identify  to  the  locations  of  the  FEP’s.  Due  to  the  difficultly  in  locating  the  energy  peaks 
in  the  spectra,  large  errors  were  associated  with  the  energy  calibrations  of  the  detectors. 
To  fix  the  problem,  the  1,274  and  511  keV  trigger  voltages  were  used  instead  of  energy 
calibrations  to  run  the  code.  The  second  problem  encountered  occurred  due  to  the  low 
activity  of  the  source.  The  oscilloscope  was  originally  set  up  so  that  it  would  only  trigger 
off  of  a  1274  keV  photon  followed  by  a  511  keV  photon.  Due  to  the  low  activity  of  the 
capillary  tube,  if  no  pulses  met  these  conditions  within  a  specified  amount  of  time,  the 
system  timed  out.  To  fix  the  problem,  the  timeout  was  increased  to  60  sec. 

Even  after  fixing  the  timeout  problem,  the  efficiency  of  the  system  was  still  poor. 
Although  a  direct  comparison  between  the  PsARS  and  PALS  results  could  no  longer  be 
made,  the  first  capillary  tube  with  the  second  capillary  tube.  The  specific  results  are 
discussed  in  the  next  section,  but  the  second  capillary  tube  did  have  a  higher  activity  than 
the  first  capillary  tube.  Also,  to  increase  the  chance  of  seeing  the  o-Ps  in  the  lifetime 
spectrum,  the  lifetime  measurement  was  conducted  with  a  bias  of  75  V  applied  to  the 
silver  layer  of  the  capillary  tube.  After  collecting  data  for  approximately  42  hours,  the 
following  results  were  obtained: 
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Lifetime  Spectrum 


Figure  40.  Lifetime  spectrum  from  measuring  the  second  capillary  tube,  with  a  bias  of  75  V  added  to 
the  silver  layer,  for  42  hours. 

As  can  be  seen  from  the  spectrum,  the  system  is  very  ineffieient  and  did  not 
colleet  enough  counts  to  run  a  PALS  fit  analysis.  However,  the  results  show  that  the 
system  is  working  properly.  As  expected,  most  of  the  counts  are  in  the  pico  second 
region  due  to  the  annihilation  of  p-Ps,  but  peaks  from  the  o-Ps  can  be  seen  starting  to 
form.  4  peaks  above  5  ns  are  forming  due  to  the  longer  lived  lifetimes  of  the  o-Ps.  These 
oecur  at  approximately  7  ns,  13  ns,  19  ns,  and  29  ns.  These  longer  lifetimes  result  from 
the  positrons  forming  o-Ps  at  the  surface  of  the  gold  capillary  tube.  Together  with  rise  in 
Sy/ly  ratio,  both  results  proved  that  at  the  capillary  tube  surface,  both  the  Sy/ly  ratio  and 
positron  lifetimes  increased.  To  achieve  a  spectrum  that  can  be  analyzed  with  PALSfit 
and  to  better  see  the  longer  lived  o-Ps,  more  counts  must  be  recorded  to  achieve  the 
required  1  million. 
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4.5  Second  Capillary  Tube 

0.3  |a,Ci  of  radioactive  Na-22  was  intended  to  be  deposited  inside  the  capillary 
tube.  However,  problems  arose  due  to  the  volume  of  solution  needed  to  deposit  the 
required  amount  of  Na-22  inside  the  capillary  tube.  When  the  first  capillary  tube  was 
created,  approximately  60  pL  of  the  radioactive  solution  was  needed  to  deposit  0.3  pCi  of 
Na-22  inside  the  tube.  However,  due  to  radioactive  decay,  this  year  85  pL  was  needed  to 
deposit  the  same  amount  of  Na-22  inside  the  capillary  tube.  The  capillary  tube  was  not 
able  to  hold  the  increased  volume  of  liquid  and  some  liquid  leaked  out  the  end  of  the 
tube.  Furthermore,  the  activity  was  concentrated  on  one  side  of  the  tube.  To  counteract 
these  problems  in  the  future,  two  modifications  may  be  made.  The  first,  the  same  Na-22 
solution  can  be  used,  but  part  of  the  water  must  be  evaporated  to  concentrate  the  Na-22. 
A  second  option  is  to  use  the  other  Na-22  solution  at  AFIT  that  has  a  higher  activity. 
After  measuring  the  second  capillary  tube  and  comparing  its  results  to  those  of  the  first 
capillary  tube,  approximately  0.424  pCi  of  Na-22  (7  Jan  2011)  was  deposited  into  the 
second  capillary  tube.  More  than  0.3  pCi  of  Na-22  was  deposited  into  the  capillary  tube 
due  to  the  precision  with  drawing  85  pL  into  the  syringe.  However,  the  result  shows  that 
heating  the  capillary  tube  on  both  sides  at  a  lower  temperature  prevents  pressure  from 
building  up  inside  that  would  have  ejected  the  solution  out. 

4.6  Future  Improvements 

Although  many  problems  with  the  PsARS  system  were  solved  from  the  previous 
year,  many  more  still  exist.  The  problem  with  the  vacuum  pressure  that  Lt  Slaughter 
encountered  was  not  resolved.  In  fact,  the  problem  got  worse.  While  Lt  Slaughter  was 
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able  to  achieve  pressures  of  10'  torr,  this  year  only  pressures  of  1.36x10’  torr  were 
achieved.  After  turning  on  the  roughing  pump,  initially  the  vacuum  pressure  dropped  to 
around  1.07x10'  torr,  but  then  over  time  the  pressure  increased  to  1.36x10’  torr.  To 
correct  the  problem,  all  of  the  metal  gaskets  in  between  the  CF  flanges  were  exchanged 
for  new  ones,  and  all  of  the  bolts  were  tightened  around  the  flanges  until  the  metal 
gaskets  could  no  longer  be  seen  in  between  the  flanges.  It  addition,  vacuum  grease  was 
added  to  the  rubber  o-rings  in  the  Klein  flanges  (KF).  For  future  experiments,  all  of 
metal  gaskets  should  be  replaced  again  and  the  o-rings  should  be  replaced  with  new 
rings.  Furthermore,  the  full  power  of  the  roughing  pump  and  sensitivity  of  the  sensor 
needs  to  be  tested  because  after  attaching  the  vacuum  sensor  directly  to  the  flange  on  the 
pump,  the  vacuum  pressure  only  reached  approximately  1x10’  torr. 

The  next  problem  encountered  was  when  the  XIA  cards  would  record  differing 
number  of  events.  Even  though  the  two  modules  would  sync  up  at  the  beginning  of  a 
run,  eventually  one  module  would  record  an  event  that  the  other  missed.  This  problem 
was  discussed  with  the  manufacturer  of  the  cards  and  they  determined  the  problem  to 
exist  within  the  DGF4C’s  firmware.  Unfortunately,  the  firmware  was  not  able  to  be 
updated,  but  the  problem  was  compensated  for  in  the  analysis  code.  In  this  experiment, 
the  problem  was  minimized  due  to  the  low  activity  of  the  source,  but  in  future 
experiments,  if  the  activity  is  increased,  the  problem  will  worsen.  The  best  solution  to  fix 
the  problem  would  be  to  buy  a  new  digitizing  card  with  more  than  four  channels  on  it.  A 
new  card  will  have  the  updated  firmware,  and  with  more  than  four  channels  on  the  card, 
the  process  of  getting  multiple  cards  to  communicate  with  each  other  is  eliminated. 
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The  next  improvement  is  to  increase  the  efficiency  of  the  analysis  code  and 
decrease  the  runtime  of  the  experiments.  In  the  analysis  code,  only  events  with  two  and 
three  hits  were  analyzed.  Events  with  three,  four,  five,  or  six  hits  could  contain  or  Sy 
photons  with  background  hits,  but  these  cases  were  ignored.  To  make  the  correction,  the 
code  would  have  to  look  for  2y  photons  in  events  with  three,  four,  five,  and  six  hits  and 
3y  photons  in  events  with  four,  five,  and  six  hits.  However,  by  looking  at  how  many 
events  in  a  measurement  contain  five  or  six  hits,  the  determination  can  be  made  to  ignore 
those  events  if  the  chance  of  seeing  one  is  sufficiently  low. 

A  fourth  improvement  can  be  made  by  correcting  the  problem  from  having  to 
delete  the  module  information  at  the  top  and  middle  of  all  of  the  Igor  Pro  output  files. 
Currently,  the  only  way  to  solve  the  problem  is  to  open  up  the  output  file  and  delete  the 
information  by  hand.  Larger  file  sizes  take  longer  to  load,  and  once  the  files  become  too 
large,  they  will  not  open.  The  problem  can  be  corrected  for  either  in  the  PsARS  3y/2y 
code  or  the  Igor  Pro  program.  First,  the  3y/2y  code  can  be  modified  to  account  for  the 
extra  information  and  pass  over  it  as  the  code  scans  the  data.  Secondly,  the  Igor  Pro 
settings  could  possibly  be  changed  using  the  GUI  to  prevent  the  output  file  from  having 
the  extra  information.  If  the  second  solution  was  chosen  and  the  activity  of  the  source 
being  measured  was  changed,  the  timing  information  on  the  modules  would  have  to  be 
checked  first  to  ensure  the  modules  are  start  up  simultaneously. 

The  last  problem  with  the  PsARS  system  occurred  due  to  the  difference  in  dead 
time  for  the  Nal(Tl)  and  BaF2  detectors.  The  Nal(Tl)  needs  to  measure  low  activity 
sources  while  the  BaF2  detectors  need  to  measure  high  activity  sources  to  prevent  the 
PALS  measurements  from  taking  a  week  to  conduct.  To  prevent  this  problem  in  the 
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future,  independent  experiments  can  be  run  with  different  activity  sources.  Although  the 
two  results  will  not  be  able  to  be  directly  compared  to  each  other,  they  will  still  provide 
useful  information  about  the  formation  and  lifetime  of  Ps. 

The  biggest  problem  encountered  with  the  PALS  system  came  from  how  the 
oscilloscope  digitized  the  data.  The  digitizing  rate  was  too  low  and  the  rounding  errors 
within  the  program  became  significant.  This  caused  the  lifetime  spectra  to  consist  on 
dirac  functions.  To  improve  the  system,  an  interpolation  technique  needs  to  be  added  to 
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the  MATLAB  PALS  code  when  searching  for  the  midrise  points  of  the  start  and  stop 
pulses.  This  will  reduce  the  rounding  error  and  improve  the  timing  resolution. 

The  other  major  problem  with  the  PALS  system  was  its  efficiency.  Besides 
measuring  a  separate,  higher  activity  source,  changes  to  the  code  need  to  be  made  to 
adjust  the  start  and  stop  modules.  With  only  one  detector  acting  as  the  start  module, 
immediately  the  efficiency  of  the  system  is  cut  in  half.  By  modifying  to  the  code  for  both 
detectors  to  act  as  the  start  modules,  the  efficiency  of  the  system  would  double. 

Future  research  on  the  capillary  tube  could  include  having  the  positron  interact 
with  different  materials.  This  can  be  done  in  one  of  two  ways.  The  first  would  be  to 
deposit  another  material  on  the  outside  of  the  gold  layer.  By  placing  a  positive  bias  on 
the  silver  layer,  the  positrons  would  be  drawn  to  the  gold  surface,  and  with  enough 
energy,  escape  from  the  gold  layer  into  new  material.  With  the  PALS  system,  lifetime 
analysis  could  be  used  to  characterize  the  size  of  the  pores  in  the  material.  A  second 
experiment  would  involve  filling  the  vacuum  chamber  with  a  gas.  Once  at  the  surface  of 
the  capillary  tube,  the  positrons  would  interact  with  the  atoms  or  molecules  in  the  gas. 
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5.  Conclusion 


A  new  3y/2y  code  was  developed  that  combined  the  advantages  of  using  both 
energy  and  geometry  filters.  These  changes  increased  the  efficiency  and  decreased  the 
runtime  of  the  system,  allowing  for  practical  measurements  to  be  conducted  on  low 
activity  sources.  The  method  was  calibrated  by  measuring  a  copper  standard  source  with 
a  known  3y/2y  ratio.  The  source  was  made  by  depositing  a  radioactive  NaCl  solution 
onto  a  copper  sample  and  evaporating  the  water.  In  previous  measurements  problems 
arose  due  to  the  capillary  bending  under  the  force  of  gravity.  This  problem  was  fixed  by 
designing  a  ring  of  polymer  that  was  placed  inside  the  vacuum  chamber  to  hold  the  free 
end  of  capillary  tube  and  keep  it  centered.  In  addition,  a  second  capillary  tube  was 
created  by  depositing  the  radioactive  NaCl  solution  within  the  tube.  The  tube  was  heated 
on  both  sides,  which  prevented  a  pressure  buildup  inside  the  tube  that  kept  the  solution 
inside. 

The  3y/2y  analysis  of  the  capillary  tube  showed  that  the  inner  electric  field 
effectively  drew  the  positrons  to  the  surface  of  the  capillary  tube.  Once  at  the  surface, 
they  formed  Ps.  With  a  bias  applied  to  the  silver  layer,  the  3yl2y  ratio  increased  from 
0.0186,  with  no  applied  field,  to  0.0208,  with  a  75  V  bias  on  the  silver  layer.  Next,  the 
bias  on  both  the  silver  and  gold  layer  was  increased  simultaneously,  but  the  voltage  drop 
between  the  two  remained  at  75  V.  The  first  electric  was  supposed  to  bring  the  positrons 
to  the  surface  of  the  capillary  tube  and  the  second  field  was  supposed  to  polarize  the  Ps, 
increasing  their  annihilation  lifetimes.  However,  the  results  showed  that  the  second 
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electric  field  did  not  polarize  the  Ps,  but  rather  forced  the  positrons  to  be  ejected  from  the 
capillary  before  they  could  form  Ps. 

After  unsuccessfully  troubleshooting  the  digital  PALS  system,  a  second  system 
was  developed,  the  MATLAB^m  PALS  system.  This  system  was  used  to  measure  the 
second  capillary  tube  with  an  inner  applied  bias  of  75  V.  Even  though  enough  counts 
were  not  collected  to  run  a  PALS  fit  analysis,  the  results  showed  that  the  new  system 
works  and  that  forcing  the  positrons  to  the  surface  of  the  gold  layer  increased  their 
lifetimes. 

This  research  lead  to  the  development  of  a  new  PALS  system  and  a  more  efficient 
PsARS  system.  Together,  the  results  from  both  systems  show  that  by  establishing  an 
electric  field  across  the  capillary  tube  wall  would  force  the  positrons  to  drop  into  the  Ore 
gap  as  they  reached  the  surface.  In  addition,  the  results  proved  that  at  the  capillary  tube 
surface  both  the  Sy/ly  ratio  and  lifetimes  increased.  Further  research  should  focus  on 
improving  the  efficiency  of  both  systems  and  either  adding  new  materials  to  the  capillary 
tube  or  filling  the  vacuum  chamber  with  different  gases  and  using  the  positrons  to 
characterize  the  materials.  This  research  provided  the  foundation  for  a  new  PALS  system 
and  led  to  the  development  of  conducting  PsARS  measurements  in  a  practical  amount  of 
time. 


83 


84 


Appendix  B.  MATLAB^m  3y/2y  Analysis  Code 


%program  determines  the  number  of  3  and  2  gamma  counts 

%save  data  as  .dat  file  in  folder,  drag  the  file  to  the  workspace  and 
%use  the  import  function  to  load  data 
data  =  A15minl250V1175V; 

[row, colm] =size (data) ; 

%initalize  variables  to  keep  track  of  the  number  of  each  type  of  count 

count2=0 ; 

count3_120=0 ; 

count3_adj  =0 ; 

count3_2=0 ; 

%Becuase  sometimes  each  module  will  end  up  will  different  event 
numbers , 

%this  will  search  for  the  smaller  event  number. 

%Find  the  number  of  events  in  the  2nd  module 
event2=data (row,  1) ; 

%find  the  number  of  events  in  1st  module 
eventl=l ; 

for  i=l : (row/2)+10 

if  eventl<data ( i ,  1 ) 
eventl=data (i, 1) ; 

end 

end 

%check  to  see  which  number  of  events  is  less 
if  eventl<event2 
event=eventl ; 

else 

event=event2 ; 

end 

%add  in  energy  calibrations 
for  i=l : event 

%calibrate  module  1  ch  1  and  2 

data  (4*i-2, 3)  =  (-2 . 04812e-7*data  (4*i-2, 3)  ^'2)  +  (5 .21453e-2*data  (4*i- 

2.3)  )  ; 

data(4*i-l,  3)  =  (2.43024e-7*data(4*i-l,3)'^2)  +  (1.98521e-2*data(4*i- 

1.3)  ) ; 

%calibrate  module  2  ch  0,1,2  and  3 

data (4*event+ ( (i-1) *4)+l,3)=(-3.31639e-8*data (4*event+ ( (i- 
1)*4)+1,3)  .  .  . 

''2)  +  (4.80379e-2*  data  (4*event+  ((i-l)*4)+l,3)); 
data (4*event+ ( (i-l)*4)+2,3)=(2.06727e-7*data (4*event+ ( (i- 
l)*4)+2,3) . . . 

^2)  +  (3.37157e-2*data (4*event+ ((i-l)*4)+2,3)); 
data (4*event+ ( (i-1) *4)+3,3)=(3.58532e-7*data (4*event+ ( (i- 
l)*4)+3,3)  . .  . 

''2)  +  (1 . 52752e-2*data  (4*event+  ((i-l)*4)+3,3)); 
data (4*event+ ( (i-1) *4)+4,3)=(-9.97343e-8*data (4*event+ ( (i- 
l)*4)+4,3)  .  .  . 

''2)  +  (3.32537e-2*data  (4*event+  ((i-l)*4)+4,3)); 
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end 


for  i=l : event 


%check;  to  see  if  all  the  times  for  the  events  match  up 
k;=0;  %use  k  to  keep  track  of  time  matches 

if  (data (4*i-3, 4) ==data (4*i-2, 4) )  &&  (data(4*i-3,4)==data(4*i- 

1,4)  )  .  .  . 


&&  (data (4*i-3 

, 4) ==data (4*i, 4) )  && 

(data (4*i-3, 4) 

==data (4*event+ ( (i-1) 

*4) 

'tl,4) ) 

&& 

(data (4*i-3, 4) 

==data (4*event+ ( (i-1) 

*4) 

't2,4) ) 

&& 

(data (4*i-3, 4) 

==data (4*event+ ( (i-1) 

*4) 

't3,4) ) 

&& 

(data (4*i-3, 4) 

==data (4*event+ ( (i-1) 

*4) 

't4,4) ) 

k=i; 

else 


for 


matching 

the 

and 


1)*4)+1,4) ) 
l)*4)+2,4) ) 
l)*4)+3,4) ) 
l)*4)+4,4) ) 


j=i-20:it20 

if  ((j>0)  &&  (j<=event2)) 

if  data(4*i-2,4) ==data ( 4*event+ ((j-l)*4)+l,4) 

%if  times  don't  match  up,  find  the  event  with  the 

%time.  check  that  all  6  times  match  because  sometimes 

%times  change  within  the  event,  even  when  the  events 

%times  match  up 

if  (data (4*i-2, 4) ==data  (4*i-l, 4) )  &&  ... 

(data (4*i-2, 4) ==data ( 4  * event t ( ( j - 

&& (data (4*i-2, 4) ==data ( 4  * event + ( ( j - 

&&  (data(4*i-2,4) ==data ( 4*event+ (  ( j - 

&&  (data(4*i-2,4) ==data ( 4*event+ ( ( j - 

k=j  ; 

end 

end 


end 


end 


end 


%if  program  could  get  the  time  from  modi  to  match  to  mod2  for  an 
event, 

%continue  on  and  look  for  coincidence 
if  k>0 

%rename  the  PMT ' s  to  make  coding  easier  to  read  and  write 
%PMT  1  is  "A",  PMT  2  is  "B" . . . 

A=data (4*i-2, 3)  ; 

B=data (4*event+ ( (k-1) *4) +1, 3) ; 

C=data (4*event+ ( (k-1) *4) +2, 3) ; 

D=data (4*i-l, 3) ; 

E=data (4*event+ ( (k-1) *4) +3, 3) ; 

F=data (4*event+ ( (k-1) *4) +4, 3) ; 
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%check  the  number  of  hits  for  each  event 
hit=0 ; 

for  j=4*i-2 : 4*i-l 
if  data ( j , 3 ) >0 
hit=hit+l ; 

end 

end 

for  j=4*event+ ( (k-1) *4) +1 : 4*event+ ( (k-1) *4) +4 
if  data ( j , 3 ) >0 
hit=hit+l ; 

end 

end 

%check  to  see  if  2  hits  are  2gamma  or  3gamma  count 
if  hit==2 
sum=0 ; 

%check  if  hits  are  on  opposite  PMT ' s  for  2gamma  count 
if  ( ( (A>400) && (D>400) )  &&  ( (A<600 ) && (D<600 ) ) ) 
sum=  A  +  D; 

elseif  (  (  (B>400) && (E>400) )  &&  ( (B<600 ) && (E<600 ) ) ) 
sum=  B  +  E; 

elseif  (  (  (0400) && (F>400)  )  &&  (  (C<600 ) && (F<600 )  )  ) 
sum3=  C  t  F; 

end 

%if  hits  are  2  511 's,  count  as  a  2  gamma  hit (resolution  at 
%  511  is  10%,  150  is  20%) 
if  ((sum>580)  &&  (sum<1124)) 
count2=count2+l ; 

end 

%check  3gamma  non-adjacent  hits  (  resolution  at  380  is  12%) 


if  ( ( (A  &&  C) >50)  && 
sum=  A  +  C; 

(  (A 

&&  C 

)<425) ) 

elseif  ( ( (B  && 
sum=  B  +  D; 

D) >50) 

&& 

(  (B 

&& 

D) <425 

elseif  ( ( (C  && 
sum=  C  t  E; 

E) >50) 

&& 

(  (C 

&& 

E) <425 

elseif  ( ( (D  && 
sum=  D  +  F; 

F) >50) 

&& 

(  (D 

&& 

F) <425 

elseif  ( ( (E  && 
sum=  E  +A; 

A) >50) 

&& 

(  (E 

&& 

A)  <425 

elseif  ( ( (F  && 

B) >50) 

&& 

(  (F 

&& 

B) <425 

sum=  F  +B; 

end 

if  ((sum>160)  &&  (sum<850))  %  resolution  at  380  at  12%, 

%  150t50  is  20% 
count3_2=count3_2+l ; 

end 

end 

%check  for  3  hits 
if  hit==3 
sum=0 ; 


87 


%check;  for  120  degree  hits  (every  other  detector)  (380  at 

12%) 

if  (  (  (A>50) && (050) && (E>50)  )  &&  ((A<425)&&  (C<425)&& 

(E<425) ) ) 

sum=  A  +  C  t  E; 

elseif ( ( (B>50 ) && (D>50 ) && (F>50 ) )  &&  ((B<425)&& 

(D<425) && (F<25) ) ) 

sum=  B  +  D  +  F; 

end 

if  ((sum>640)  &&  (sum<1124))  %  340+340  at  12%  +50  at  20% 
count3_120=count3_120+l ; 

end 

%check  for  opposing  hits  with  an  adajecent  hit,  380  at  12%) 
if  (  (((A  Sc  Sc  D)&&(  C  I  E  )  )>50)  &&  ((A  &&  C  &&  D  && 

E)  <425) ) 

sum=  A+E+D+C; 

elseif  (  ( ( (B  &&  E) && (  D  |  F  )  )>50)&&((B  &&  E  &&  D  && 

F)  <425) ) 

s  um=  B  +  D  +  E  +  F; 

elseif  (  ( ( (C  &&  F) && (  E  |  A  )  )>50)&&((A  &&  C  &&  F  && 

E) <425) ) 

sum=  C  +  E  +  F  +  A; 

elseif  (  ( ( (D  &&  A) && (  F  |  B  )  )>50)&&((A  &&  F  &&  D  && 

B) <425) ) 

s  um=  D  +  F  +  A  +  B; 

elseif  (  (((E  &&  B)&&(  A  |  C  )  )^50)&&((A  &&  C  &&  B  && 

E) <425) ) 

sum=  E+A+B+C; 

elseif  (  (((F  &&  C)&&(  B  |  D  )  )>50)&&((F  &&  C  &&  D  && 

B) <425) ) 

s  um=  F  +  B  +  C  +  D; 

end 

if  ((sum>308)  &&  (sum<1124))  %  150+150+50  at  12% 
count3_adj  =count3_adj  +1 ; 

end 

end 

end 

end 

count2 

count3Tot=count3_12  0+count3_adj  +count3_2 
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Appendix  C.  data_cal_stream.m  File 

%  Create  a  TCPIP  interface  and  device  object  -  Works  with  4  Channels 
clc 

clear  all 
close  all 

Run_Start_time  =  fix (clock)  %  get  time  for  starting  acquisition 
dpo . InputBuf ferSize  =  300000;  %buffer  size 

dpo  =  instrf ind ( ' Type '  ,  'tcpip',  ' RemoteHost ' ,  '169.254.235.51', 

' RemotePort ' ,  4000,  'Tag',  ' ' ) ;  %Identifies  object  connecting  to 

if  isempty(dpo) 

dpo  =  tcpip (' 169 . 254 . 235 . 51 '  ,  4000);  %  ip  address  of  DPO7104 

else 

fclose (dpo) ; 
dpo  =  dpo  ( 1 ) ; 

end 


deviceObj  =  icdevice ( ' tektronix_dpo7104_rl .mdd ' ,  dpo) 
scope  using  driver  called  tektronix  dpo  7104  rl.mdd 


%Connecting  to 


%Connect  to  the  scope 
connect (deviceObj ) ; 

%%  Scope  Settings 
fprintf (dpo, ' Select : Ch3  On') 
fprintf (dpo, ' Select : Ch4  On') 
fprintf (dpo, ' CH3 : Position  0' 
fprintf (dpo, ' CH4 : Position  O' 


turns  input  for  channel  3 
turns  input  for  channel  4 
5  Aligns  channel  3  at  0  V 
5  Aligns  channel  4  at  0  V 


on 

on 


Places  signals  at  designated 
Sets  vertical  scale  for  channel  3 
Sets  vertical  scale  for  channel  4 


Toggles  enhanced  triggering  mode 
Sets  the  A  trigger  mode  as  a 


fprintf (dpo, ' Horizontal : Position  20') 
horizontal  location  (10) 
fprintf (dpo,  ' ch3 : Scale  .  04  '  ) ; 

(V/div) 

fprintf (dpo, ' ch4 : Scale  . 04 ' ) ; 

(V/div) 

%  Trigger  Setting 

fprintf (dpo,  ' Trigger : Enhance  Off ' )  ; 
fprintf (dpo, ' Trigger : A: Type  Logic ' ) , 

Logic 

fprintf (dpo, ' Trigger : A : Logic : Class  Pattern ' ) ; 
use  the  Pattern  Mode 

fprintf (dpo,  ' Trigger : A: Logic : Function  And ' ) ; % 
boolean  logic  to  perform  on  triggered  events 
%  AND  specifies  to  trigger  if  all  conditions  are  true. 

%  NANd  specifies  to  trigger  if  any  of  the  conditions  are  false. 

%  NOR  specifies  to  trigger  if  all  conditions  are  false. 

%  OR  specifies  to  trigger  if  any  of  the  conditions  are 
fprintf (dpo,  ' Trigger : A: Logic : Pattern : Input : CH3  High ' ) ; 

Sets  importance  level  for  channel  1  (High) 
fprintf (dpo, ' Trigger : A: Logic : Pattern : Input : CH4  High ' ) ; 
importance  ievel  for  channel  2  (High) 
fprintf (dpo,  ' Trigger : A: Logic : Input : Format  Binary ' ) ; 
to  Binary  Can  also  be  'Hexadecimal' 

fprintf (dpo,  ' Trigger : A: Logic : Pattern : When  False');  % 

{ TRUe I FALSe I LESSThan 1 MOREThan }  Sets  condition  for  meeting  boolean  logic 


%  Sets  the  A  trigger  to 
{AND  I NANd I  NOR  I  OR}  Sets 


true . 

‘  %  (HIGH  I  LOW  I  X} 

’  %  Sets 

Sets  data  format 
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fprintf (dpo, ' Trigger : A : Logic : Threshold : CH3  -.01');  %  Sets  trigger  level 
for  Channel  3 

fprintf (dpo, ' Trigger : A: Logic : Threshold : CH4  -.01');  %  Sets  trigger  level 
for  Channel  4 

fprintf (dpo, ' Trigger :A:Mode  Normal');  %  Sets  trigger  mode  to  Normal 
(Other  option  is  Auto  ->  forces  trigger  at  timeout  period) 
fprintf (dpo,  ' ACQUIRE : MODE  HiRes ' )  ;  %Sets  sampling  mode  to  High 
Resolution  Mode 


NumberOfPoints  =  str2num (query (dpo,  ' HORIZONTAL : MODE : RECORDLENGTH? ')) ;  % 
Queries  the  record  length  (number  of  points  per  waveform) 

SampleRate  =  str2num (query (dpo, ' HORIZONTAL : MODE : SAMPLERATE? ')) ;  % 
Queries  the  sample  rate  of  the  scope 

TimeSpan  =  NumberOf Points/SampleRate;  %  Returns  seconds  per  point 
timearray  =  linspace ( 0 , TimeSpan, NumberOf Points ) ;  %  Array  of  time 


fprintf (dpo,  ' Data ; Encdg  Ribinary');  %  Specifies  signed  integer-data 
point  representation  with  the  most  significant  byte  transferred  first 
(Try  SIRBINARY) 

fprintf (dpo, ' WFMOutpre : BYT_NR  2');  %  Sets  the  binary  width  of  the  data 
fprintf (dpo,  'DATA: START  200');  %  Sets  the  beginning  of  the  waveform 
data  to  transfer 

fprintf (dpo, 'DATA: STOP  1450');  %  Sets  the  end  of  the  waveform  data  to 
transfer 

fprintf (dpo, ' WFMOutpre : ENCDG  BIN')  %  Sets  incoming  data  to  binary 
format 

fprintf (dpo, ' WFMInpre : ENCDG  BINARY')  %  Sets  outgoing  data  to  binary 
format  (Lose?) 

fprintf (dpo, ' WFMInpre : BN_FMT  RI ' ) ;  %  Specifiies  integer-data  point 
representaoitn  (Lose?) 

fprintf (dpo, ' WFMInpre : BYT_NR  2');%set  2  byte  format  =  intI6  (Lose?) 
fprintf (dpo, ' Horizontal : FastFrame : State  O'); 

fprintf (dpo, 'DATA: SOURCE  Ch3 ' ) ;  %  Sets  data  source  to  Channel  3 
Ch3Ymult  =  str2num (query (dpo, ' WFMOutpre : YMult? ')) ;  %  Returns  vertical 

scale  factor 

Ch3Yoff  =  str2num (query (dpo, ' WFMOutpre : Yzero? ')) ;  %  Returns  the  offset 

factor 

fprintf (dpo, 'DATA: SOURCE  Ch4 ' ) ;  %  Sets  data  source  to  Channel  4 
Ch4Ymult  =  str2num (query (dpo, ' WFMOutpre : YMult? ')) ;  %  Returns  vertical 
scale  factor 

Ch4Yoff  =  str2num (query (dpo, ' WFMOutpre : Yzero? ')) ;  %  Returns  the  offset 

factor 


fprintf (dpo, 'DATA: SOURCE  Ch3,CH4');  %  Reinitiates  data  source  to  all 

channels 

nSamples=5500000 ;  %  Sets  the  number  of  sammples  to  acquire  (set  really 
large  to  allow  for  time  acquisition  instead  of  number  of  acquisitions) 
minutes_run  =1;  %  Sets  how  many  minutes  to  acuire  data  for 

seconds_run  =  minutes  run*60;  %  Changes  length  of  acquisition  to 

seconds 

tstart  =  tic;  %  Start  acquistion  timer 
for  i  =  l:nSamples 

fprintf (dpo, 'CURVE?');  %  Call  waveform  (writes  out  in  list  mode 

Channel  1  first  then  Channel?  then...) 


90 


Channel  (:,1)  =  Ch3Ymult*binblock;read  (dpo,  '  inti  6  '  ) +Ch3Yof  f ;  %  Reads 
Channel  3  and  scales  it 

Channel (:,  2)  =  Ch4Ymult*binblockread (dpo,  ' inti  6 ' ) +Ch4Yof f ;  %  Reads 
Channel  4  and  scales  it 

fprintf  (dpo,  '  *wai  '  )  %  Makes  scope  wait  to  finish  tranfer  before 

executing  next  command 

fread (dpo, 1 ) ;  %  Clears  extra  terminator  from  'CURVE?'  cmd 
fprintf (dpo, ' *wai ' ) ;  %  Makes  scope  wait  to  finish  tranfer  before 
executing  next  command 
if  i==l 

save ( ' Na22_Spec2 . txt '  ,  ' Channel ',  '-ascii',  '-tabs') 

else 

save ( ' Na22  Spec2 . txt ' , ' Channel ' , ' -append ' , ' -ascii ' , ' -tabs ' ) 

end 

if  toc(tstart)  >  seconds_run  %  Exit  for  loop  when  timing  criteria 
has  been  met 
break 

end 

end 

Sample_Collection_Time  =  toc(tstart)  %  Gives  how  long  acquisition  took 
Run_Stop_time  =  fix (clock)  %  Gives  current  time 


91 


Appendix  D.  cal_scan.m  File 


close  all 
clear  all 
clc 

Run_Start_time  =  fix (clock) 
clear  segarray; 

NumberOf Points  =  1251; 

block_size  =  NumberOfPoints*100;  %must  match  the  number  of  points  per 
waveform  you  transmitt 

%  describe  the  format  of  the  data 

%  for  more  information,  see  the  textscan  reference  page 

file_id  =  f open ( ' Na22_Spec . txt ' ) ;  %establish  ID  for  file  you  open  up 

count  =  0; 
channel3=  [  ] ; 
channel4  =  [  ] ; 

while  ~f eof ( f ile_id)  %runs  to  end  of  file 
count  =  count+1; 

segarray  =  textscan (file_id, ' %n%n block_size, ' delimiter ; 
%reads  in  a  block  of  data 

[channel3_new, channel4_new]  =  cal_func ( segarray, NumberOf Points ) ; 
channel3  =  [channel3;  channel3_new] ; 
channel4  =  [channel4;  channel4  new]; 

end 

fclose(file  id) ; 


[channel3_ph, xoutch3] =hist (channel3 ( :  ,  1)  ,  1024)  ; 
[channel4_ph, xoutch4 ] =hist (channel 4 (:,1),1024); 

figure  ( 1 ) 

plot (xoutch3, channel3_ph) 
title ( ' ch3 ' ) 

figure  (2 ) 

plot (xoutch4, channel4_ph) 
title ( ' ch4 ' ) 

Run_Stop_time  =  fix (clock) 
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Appendix  E.  cal_func.m  Function  File 


function  [channel3_out, channel4_out]  = 
cal_func ( segarray, NumberOf Points ) 

Channels  =  cell2mat ( segarray (:, 1 )) ; 

Channels  =  reshape (Channels , NumberOf Points ,[]) ; 

Channel4  =  cell2mat ( segarray (:,  2 ))  ; 

Channel4  =  reshape (Channel4 , NumberOf Points , [ ] ) ; 

clear  segarray 

SampleRate  =  1 . OOOOe+009;  %  Sample  rate  of  scope 

TimeSpan  =  NumberOf Points/SampleRate;  %  Returns  seconds  per  point 
timearray  =  linspace ( 0 , TimeSpan, NumberOf Points ) ;  %  Array  of  time 
detnum=2;  %  Number  of  detectors 

trigger_level  =  lOe-S;  %  Trigger  level  in  V,  must  be  above  noise 
[num_points , num_wfms ]  =  size (Channels) ;  %  Determines  how  many  waveforms 
were  recorded 

ChannelSfilt  =  sgolayfilt (Channels, 9, IS) ; 

Channel4filt  =  sgolayfilt (Channel4, 9,  IS)  ; 

clear  Channels 
clear  Channel4 

[pulse_peak_chS , min_index_chS ]  =  min (ChannelSfilt) ;  %  Find  the  pulse 

height  of  each  pulse  and  the  index  of  that  pulse  height  for  channel  S 
[pulse_peak;_ch4 ,  min_index_ch4  ]  =  min  (Channel4f  ilt)  ;  %  Find  the  pulse 

height  of  each  pulse  and  the  index  of  that  pulse  height  for  channel  4 

midrise  time  chS  =  zeros (l,num  wfms); 
midrise_time_ch4  =  zeros ( 1 , num_wfms ) ; 
pulse_order_time_comp  =  zeros (num_wfms , 2 ) ; 
pulse_order_index  =  zeros (num  wfms, 2); 
energy_order_index  =  zeros (num_wfms , 2 ) ; 

%  for  i  =  1 : num_wfms 

q, 

o 

[midrise_time_chl (i) ] =f ind_midrise_time (Channellf ilt ( : , i) , timearray, tri 
gger_level ) ;  %  Calls  function  " f ind_trigger_time"  and  gives  the  timming 
for  each  triggered  event  for  Channel  S  (time  is  off  of  the  trigger) 

q, 

0 

[midrise_time_ch2 (i) ] =f ind_midrise_time (Channel2f ilt ( : , i) , timearray, tri 
gger_level ) ;  %  Calls  function  " f ind_trigger_time"  and  gives  the  timming 
for  each  triggered  event  for  Channel  4  (time  is  off  of  the  trigger) 

%  end 

pulse_peak_ch3  =  -pulse  peak_ch3;  %  Changes  pulse  values  from  negative 
to  positive  for  channel  3 
en_cal_ch3  =  pulse  peak_ch3; 

pulse_peak_ch4  =  -pulse_peak_ch4 ;  %  Changes  pulse  values  from  negative 

to  positive  for  channel  4 
en_cal_ch4  =  pulse  peak_ch4; 
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for  j  =  l:num_wfms 

if  en_cal_ch3 ( j )  <  trigger_level  %  If  waveform 
trigger  level  then  call  it  a  NaN 

en_cal_ch3 ( j )  =  NaN;  %  If  waveform  peak  is 
then  call  it  a  NaN 
end 

if  en_cal_ch4 ( j )  <  trigger_level  %  If  waveform 
trigger  level  then  call  it  a  NaN 

en_cal_ch4 ( j )  =  NaN;  %  If  waveform  peak  is 
then  call  it  a  NaN 
end 


peak  is  below 
below  trigger 

peak  is  below 
below  trigger 


level 


level 


end 

[ch3row, ch3col]  =  size (en_cal_ch3) ; 
[ch4row, ch4col]  =  size (en_cal_ch4 ) ; 
if  ch3col>0 

channel3_out  =  en__cal_ch3  '  ; 

else 

channel3_out  =  []; 

end 


if  ch4col>0 

channel4_out  =  en_cal_ch4 ' ; 

else 

channel4_out  =  []; 

end 


end 
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Appendix  F.  Stream_data_PALS.m  File 

%  Create  a  TCPIP  interface  and  device  object  -  Works  with  4  Channels 
clc 

clear  all 
close  all 

Run_Start_time  =  fix (clock)  %  get  time  for  starting  acquisition 
dpo . InputBuf ferSize  =  3000000;  %buffer  size 

dpo  =  instrf ind ( ' Type '  ,  'tcpip',  ' RemoteHost ' ,  '169.254.235.51', 

' RemotePort ' ,  4000,  'Tag',  ' ' ) ;  %Identifies  object  connecting  to 

if  isempty(dpo) 

dpo  =  tcpip (' 169 . 254 . 235 . 51 '  ,  4000);  %  ip  address  of  DPO7104 

else 

f close (dpo) ; 
dpo  =  dpo  ( 1 ) ; 

end 


deviceObj  =  icdevice ( ' tektronix_dpo7104_rl .mdd ' ,  dpo) 
scope  using  driver  called  tektronix_dpo_7104_rl .mdd 


%Connecting  to 


%Connect  to  the  scope 
connect (deviceObj ) ; 
dpo. Timeout  =  60; 

%%  Scope  Settings 
fprintf (dpo,  ' Select : Ch3  On') 
fprintf (dpo, ' Select : Ch4  On') 
fprintf (dpo, ' CH3 : Position  O' 
fprintf (dpo, ' CH4 : Position  O' 


turns  input  for 
turns  input  for 
^  Aligns  channel 
5  Aligns  channel 


fprintf (dpo, ' Horizontal : Position  20') 
horizontal  location  (10) 
fprintf (dpo, ' ch3 : Scale  .05'); 

(V/div) 

fprintf (dpo, ' ch4 : Scale  . 05 ' ) ; 

(V/div) 

%  Trigger  Setting 

fprintf (dpo,  ' Trigger : Enhance  Of f ' )  ; 

'  Trigger : A: Type  Logic') 


channel  3  on 
channel  4  on 

3  at  0  V 

4  at  0  V 

Places  signals  at  designated 


Sets  vertical  scale  for  channel  3 


Sets  vertical  scale  for  channel  4 


Toggles  enhanced  triggering  mode 
Sets  the  A  trigger  mode  as  a 


fprintf (dpo. 

Logic 

fprintf (dpo,  ' TRIGger : A: LOGIC : SETHold : CLOCk : EDGE  FALL' ) ; 

fprintf (dpo,  ' TRIGger : A : LOGIC : SETHold : CLOCk : SOUrce  CH3 ' ) ; 

fprintf (dpo, ' TRIGGER: A: LOGIC : SETHOLD : CLOCK ; THRESHOLD : Ch3  -0.06' ) ; 

fprintf (dpo,  ' TRIGGER: A: LOGIC : SETHOLD : DATA: SOURCE  CH4 ' )  ; 

fprintf (dpo,  ' TRIGGER: A: LOGIC : SETHOLD : DATA: THRESHOLD : CH4  -0 . 01  '  )  ; 

' TRIGGER:A: LOGIC: SETHOLD: HOLDTIME  0.00'); 
'TRIGger:A:LOGIc:SETHold:QUAlify  OCCurs ' ) ; 
'TRIGGER:A:LOGIC:SETHOLD:SETTIME  IOOE-9' ) ; 

%  fprintf (dpo, ' Trigger : A: Logic : Class  Pattern');  %  Sets  the  A  trigger 
to  use  the  Pattern  Mode 

%  fprintf (dpo,  ' Trigger : A: Logic : Function  Or') 

%  fprintf (dpo, ' Trigger : A: Logic : Pattern : Input : CH3  High') 

Sets  importance  level  for  channel  3  (High) 

%  fprintf (dpo, ' Trigger : A: Logic : Pattern : Input : CH4  High') 
importance  level  for  channel  4  (High) 


fprintf (dpo, 
fprintf (dpo, 
fprintf (dpo. 


{HIGH  I  LOW  I  X} 


Sets 
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Sets  data  format 


fprintf (dpo,  ' Trigger : A: Logic : Input : Format  Binary');  % 
to  Binary  Can  also  be  'Hexadecimal' 

%  fprintf (dpo, ' Trigger : A: Logic : Pattern : When  True');  % 

{ TRUe I FALSe I LESSThan I MOREThan }  Sets  condition  for  meeting  boolean  logic 
%  fprintf (dpo,  ' Trigger : A: Logic : Threshold : CH3  -.095');  %  Sets  trigger 
level  for  Channel  3 

%  fprintf (dpo, ' Trigger : A: Logic : Threshold : CH4  -.035');  %  Sets  trigger 
level  for  Channel  4 

fprintf (dpo, ' Trigger :A:Mode  Normal');  %  Sets  trigger  mode  to  Normal 
(Other  option  is  Auto  ->  forces  trigger  at  timeout  period) 
fprintf (dpo, 'ACQUIRE :MODE  HiRes ' ) ;  %Sets  sampling  mode  to  High 
Resolution  Mode 


NumberOfPoints  =  str2num (query (dpo, ' HORIZONTAL : MODE : RECORDLENGTH? ')) ;  % 
Queries  the  record  length  (number  of  points  per  waveform) 

SampleRate  =  str2num (query (dpo, ' HORIZONTAL : MODE : SAMPLERATE? ') )  % 

Queries  the  sample  rate  of  the  scope 

TimeSpan  =  NumberOf Points/SampleRate;  %  Returns  seconds  per  point 
timearray  =  linspace ( 0 , TimeSpan, NumberOf Points ) ;  %  Array  of  time 


fprintf (dpo, ' Data : Encdg  Ribinary');  %  Specifies  signed  integer-data 
point  representation  with  the  most  significant  byte  transferred  first 
(Try  SIRBINARY) 

%  fprintf (dpo, ' Data:Width  2'); 

fprintf (dpo, ' WFMOutpre : BYT_NR  2');  %  Sets  the  binary  width  of  the  data 
fprintf (dpo, 'DATA: START  200');  %  Sets  the  beginning  of  the  waveform 
data  to  transfer 

fprintf (dpo, 'DATA: STOP  700');  %  Sets  the  end  of  the  waveform  data  to 
transfer 

fprintf (dpo, ' WFMOutpre : ENCDG  BIN')  %  Sets  incoming  data  to  binary 
format 

fprintf (dpo, ' WFMInpre : ENCDG  BINARY')  %  Sets  outgoing  data  to  binary 
format  (Lose?) 

fprintf (dpo, ' WFMInpre : BN_FMT  RI ' ) ;  %  Specifiies  integer-data  point 
representaoitn  (Lose?) 

fprintf (dpo, ' WFMInpre : BYT_NR  2');%set  2  byte  format  =  intIG  (Lose?) 
fprintf (dpo, ' Horizontal : FastFrame : State  O'); 

fprintf (dpo,  'DATA: SOURCE  Ch3 ' )  ;  %  Sets  data  source  to  Channel  3 
Ch3Ymult  =  str2num (query (dpo, ' WFMOutpre : YMult? ')) ;  %  Returns  vertical 

scale  factor 

Ch3Yoff  =  str2num (query (dpo, ' WFMOutpre : Yzero? ')) ;  %  Returns  the  offset 

factor 

fprintf (dpo, 'DATA: SOURCE  Ch4 ' ) ;  %  Sets  data  source  to  Channel  4 
Ch4Ymult  =  str2num (query (dpo, ' WFMOutpre : YMult? ')) ;  %  Returns  vertical 
scale  factor 

Ch4Yoff  =  str2num (query (dpo, ' WFMOutpre : Yzero? ')) ;  %  Returns  the  offset 

factor 


fprintf (dpo, 'DATA: SOURCE  Ch3,CH4');  %  Reinitiates  data  source  to  all 

channels 

nSamples=1000000000 ;  %  Sets  the  number  of  samples  to  acquire  (set 
really  large  to  allow  for  time  acquisition  instead  of  number  of 
acquisitions ) 

minutes_run  =  1440;  %  Sets  how  many  minutes  to  acquire  data  for 
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seconds_run  =  minutes  run*60;  %  Changes  length  of  acuisition  t 

o  seconds 

tstart  =  tic;  %  Start  acquistion  timer 
for  i  =  IrnSamples 

fprintf (dpo, 'CURVE?');  %  Call  waveform  (writes  out  in  list  mode 

Channel  1  first  then  Channel2  then...) 

Channel (:,1)  =  Ch3Ymult*binblockread (dpo,  ' inti 6 ' ) +Ch3Yof f ;  %  Reads 
Channel  3  and  scales  it 

Channel (:,  2)  =  Ch4Ymult*binblockread (dpo,  ' inti  6 ' ) +Ch4Yof f ;  %  Reads 
Channel  4  and  scales  it 

fprintf (dpo,  ' *wai ' ) ;  %  Makes  scope  wait  to  finish  tranfer  before 
executing  next  command 

fread (dpo, 1 ) ;  %  Clears  extra  terminator  from  'CURVE?'  cmd 
fprintf (dpo, ' *wai ' ) ;  %  Makes  scope  wait  to  finish  tranfer  before 
executing  next  command 

pulse_peak_ch3  =  min (Channel (:, 1 )) ;  %  Find  the  pulse  height  of 

each  pulse  for  channel  3 

pulse_peak_ch4  =  min (Channel (:, 2 )) ; 
if  (pulse_peak_ch3>- . 2  &&  pulse_peak_ch4>- . 2 ) 
if  i==l 

save ( ' C7  5V_Cap2 . txt ' ,  ' Channel ',  '-ascii',  '-tabs') 

else 

save ( ' C75V_Cap2 . txt ' , ' Channel ' , ' -append ', '-ascii', '-tabs') 

end 

end 

if  toe (tstart)  >  seconds_run  %  Exit  for  loop  when  timing  criteria 
has  been  met 
break 

end 

end 

Sample_Collection_Time  =  toe (tstart)  %  Gives  how  long  acquisition  took 
Run_Stop_time  =  fix (clock)  %  Gives  current  time 
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Appendix  G.  PALS_meas.m  File 


close  all 
clear  all 
clc 

Run_Start_time  =  fix (clock) 
clear  segarray; 

NumberOf Points  =  501;  %must  match  the  number  of  points  per  waveform  you 
transmitt 

block_size  =  NumberOfPoints*!; 

%  describe  the  format  of  the  data 

%  for  more  information,  see  the  textscan  reference  page 

file_id  =  f open ( ' B75V_Cap2 . txt ' ) ;  %establish  ID  for  file  you  open  up 

count  =  0; 

PALS_time=  [ ] ; 

while  ~f eof ( f ile_id)  %runs  to  end  of  file 
count  =  counttl; 

segarray  =  textscan (file_id, ' %n%n block_size, ' delimiter ; 

%reads  in  a  block  of  data 

[ PALS_time_new]  =  pals_func ( segarray, NumberOf Points ) ; 

PALS_time  =  [PALS_time;  PALS_time_new] ; 

end 

fclose(file  id) ; 


[pals_time, xoutPals] =hist (PALS_time ( :  ,  1)  ,  1024)  ; 
figure  ( ) 

plot (xoutPals, pals_time) 
title (' Lifetime  Spectrum') 
xlabel('time  (sec) ') 
ylabel ( ' Counts ' ) 

Run_Stop_time  =  fix (clock) 
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Appendix  H:  PALS_func.m  File 


function  [PALS_time]  =  pals_func ( segarray, NumberOf Points ) 

PALS_time=[] ; 

Channels  =  cell2mat ( segarray (:,  1 ))  ; 

Channels  =  reshape (Channels , NumberOf Points ,[])  ; 

Channel4  =  cell2mat ( segarray (:, 2 )) ; 

Channel4  =  reshape (Channel4 , NumberOf Points , [ ] ) ; 

clear  segarray 

SampleRate  =  l.OOOOe+10;  %  Sample  rate  of  scope 

TimeSpan  =NumberOf Points/SampleRate;  %  Returns  seconds  per  point 
timearray  =  linspace ( 0 , TimeSpan, NumberOf Points ) ;  %  Array  of  time 
detnum=2;  %  Number  of  detectors 

max_coincidence_window  =  lOOe-9;  %  Width  of  coincidence  window  in 
seconds 

trigger_level  =  lOe-S;  %  Trigger  level  in  V,  must  be  above  noise 
[num_points , num_wfms ]  =  size (Channels) ;  %  Determines  how  many  waveforms 
were  recorded 

ChannelSfilt  =  sgolayfilt (Channels,  9,  S9)  ; 

Channel4filt  =  sgolayfilt (Channel4, 9, S9) ; 

[pulse_peak_chS , min_index_chS ]  =  min (ChannelSfilt) ;  %  Find  the  pulse 

height  of  each  pulse  and  the  index  of  that  pulse  height  for  channel  S 
[pulse_peak;_ch4 ,  min_index_ch4  ]  =  min  (Channel4f  ilt)  ;  %  Find  the  pulse 

height  of  each  pulse  and  the  index  of  that  pulse  height  for  channel  4 

trigger  level_chS  =  .01;  %V 
trigger_level_ch4  =  .01;  %V 

midrise_time_chS  =  zeros ( 1 , num_wfms ) ; 
midrise_time_ch4  =  zeros ( 1 , num_wfms ) ; 
pulse  order_time_comp  =  zeros (num_wfms , 2 ) ; 
pulse_order_index  =  zeros (num_wfms , 2 ) ; 
energy_order_index  =  zeros (num_wfms , 2 ) ; 

for  i  =  l:num_wfms 

[midrise_time_chS ( i ) ] =f ind_midrise_time (ChannelSfilt ( : , i ) , timearray, tri 
gger_level);  %  Calls  function  " f ind_trigger_time"  and  gives  the  timming 
for  each  triggered  event  for  Channel  S  (time  is  off  of  the  trigger) 

[midrise_time_ch4 (i) ] =f ind_midrise_time (Channel4f ilt ( : , i) , timearray, tri 
gger_level);  %  Calls  function  " f ind_trigger_time"  and  gives  the  timming 
for  each  triggered  event  for  Channel  4  (time  is  off  of  the  trigger) 
end 

pulse_peak_chS  =  -pulse_peak_chS ;  %  Changes  pulse  values  from  negative 

to  positive  for  channel  S 
en_cal_chS=pulse_peak_chS ; 
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pulse_peak;_ch4  =  -pulse  peak;_ch4;  %  Changes  pulse  values  from  negative 
to  positive  for  channel  4 
en_cal_ch4=pulse_peak_ch4 ; 

for  j  =  l:num_wfms 

if  en_cal_ch3 ( j )  <  trigger_level_ch3  %  If  waveform  peak  is  below 
trigger  level  then  call  it  a  NaN 

en_cal_ch3 ( j )  =  NaN;  %  If  waveform  peak  is  below  trigger  level 
then  call  it  a  NaN 

midrise_time_ch3 ( j )  =  NaN;  %  If  waveform  peak  is  below  trigger 
level  then  call  it  a  NaN 
end 

if  en_cal_ch4 ( j )  <  trigger_level_ch4  %  If  waveform  peak  is  below 
trigger  level  then  call  it  a  NaN 

en_cal_ch4 ( j )  =  NaN;  %  If  waveform  peak  is  below  trigger  level 
then  call  it  a  NaN 

midrise_time_ch4 ( j )  =  NaN;  %  If  waveform  peak  is  below  trigger 
level  then  call  it  a  NaN 
end 


Duplicate 

Duplicate 


pulse_order_time_comp ( j , 1 )  =  midrise_time_ch3 ( j ) 
timing  values  to  place  in  array  to  sort  pulses 

pulse_order_time  comp (j, 2)  =  midrise_time_ch4 ( j ) 
timing  values  to  place  in  array  to  sort  pulses 
lowerbound_511  =  .01; 
upperbound_511  =  .056; 
lowerbound_1274  =  .06; 
upperbound_1274  =  .15; 
if  (  (midrise_time_ch4 ( j ) - 
midrise_time_ch3 ( j ) <max_coincidence_window) 
midrise_time  ch3(j)>0)  ) 

if  (en_cal_ch4<upperbound_511  &&  en_cal_ch4>lowerbound_511)  && 
(en_cal_ch3<upperbound_1274  &&  en_cal_ch3>lowerbound_1274) 

PALS_time(j)  =  abs (midrise_time_ch3 ( j ) - 
midrise_time_ch4 ( j ) ) ; 

else 

%  PALS_time(j)  =  []; 

end 

end 

end 


&&  (midrise  time  ch4(j) 


end 
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Appendix  I.  find_midrise_time.m  File 


%%Find  the  time  for  halfway  up  rise  of  pulse 
function  [halfpeak_time]  = 

f ind_midrise_time (Channel, timearray, trigger_level ) 
if  min (Channel)  <=  -trigger_level 
percentSO  =  . 5*min (Channel ) ; 
midrise  =  percentSO; 

indexabove  =  find (Channel<midrise, 1 first ' ) 
if  indexabove  >  1 

indexbelow  =  indexabove  -  1; 

else 

indexbelow  =  1; 

end 

abovediff  =  abs (Channel (indexabove) -midrise) 
belowdiff  =  abs (Channel (indexbelow) -midrise) 
if  abovediff  <  belowdiff 

halfpeak_time  =  timearray (indexabove) ; 

else 

halfpeak_time  =  timearray ( indexbelow) ; 

end 

else 

halfpeak_time  =  NaN; 

end 

end 
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